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ON THE SURFACE TENSION OF ROTATING LIQUIDS* 
By M. BORNEAS AND I. BABUTIA 


Physical Laboratory; Polytechnic Institute, Timisoara, Roumania 
( Received June 10, 1960 revised manuscript received October 23, 1960) 


The surface tension of liquids changes if they are put in rotation. This is an effect specific 
to rotation — a rotational kinetic effect (r. е.). The principal factors determining the г. e. are: 
the composition of the liquid, its temperature, the velocity of rotation, and its previous treat- 
ment. The experimental procedure is described. The dependence of the r. e. on the tempera- 
ture is shown, and values for different frequencies are given. The previous thermal and mechan- 
ic treatment is discussed. It is shown that the т. e. vanishes at a certain temperature characte- 
ristic of the liquid, but independent of the velocity of rotation. Finally, an attempt at explaining 


the phenomenon is outlined. 


1. Introduction 


Measuring with a Lecomte du Noüy ring tensiometer, one of us [1] showed that the 
surface tension of rotating liquids at low temperatures presents values differing from those 
of liquids at rest. It has been pointed out [2] that this phenomenon is not identical with the 
effect described by Piccardi et al. [3], but is specific to rotation, namely, it is a rotational 
kinetic effect (r. e.). It has also been shown [1] that the r. e. does not depend on the relative 
rotation of the liquid with respect to the ring, but on the rotation with respect to the labora- 
tory. 

The principal factors determining the r. e. are: the composition of the liquid, its 
temperature, its rotation velocity, and previous treatment (mechanical and thermal). 

In order to investigate the influence of these factors on the r. e., we appropriately adapt- 
ed the Lecomte du Noüy tensiometer. We used a torsion wire adequately prepared [4] 
The dish with the liquid investigated was surrounded with a double-walled vessel containing 
water or xylene. The aim of this vessel was, on the one hand, to heat the liquid, and on 
the other, to keep the, temperature constant for a longer period of time. Cooling of the 
liquid was achieved by introducing cold water or xylene into the vessel, and for low 
temperatures, by placing the liquid in a refrigerator. 

We measured the temperature with a thermometer placed at the centre of the ring. 
This position was made possible by constructing the suspension fork-wise. With this arrange- 


ж This paper is a synthetical report of experimental results previously published, together with new results. 
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ment, the results were not affected by the presence of the thermometer, as shown by 
previous experiments [5]. 

Maximum error of the breaking angle for different liquids ranged from 0,09% to 
0.22%, which corresponds to an error of the value of the surface tension between 0.14% 
and 0.27%. 


П. Temperature dependence 


Measurements of the surface tension of rotating liquids at different temperatures[2], 
[6], [7]. [8]. [9]. [10]. [11]. led to the conclusion that the temperature-dependence curves, 
from freezing point up to boiling point, have the general form of Fig. 1. 
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Fig. 1. The general form of the curves showing the temperature dependence of the rotational-kinetic surface 
tension. 
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Fig. 2. The temperature dependence of the static (@) and rotational-kinetic (О) surface tension for butyl alcohol 
between 9 and 19°С. 
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Parts a-b and c-e of the curve are situated above the straight line given in the literature 
for liquids at rest [12]. Part b-c exhibits a minimum, which in some cases remains above 
_ the rest line (as е. g. for butyl alcohol, Fig. 2), and in other cases passes beneath, i. e. the 
r. e. becomes negative (as e. g. for glycol, Fig. 7). 

Part d-e shows a decrease of the r. e., and part e-f coincides with tbe rest line, namely, 
in the high temperature range the r. e. appears no longer. 

The curves of the foregoing type (Fig. 1) are reproducible in identical conditions. 
These conditions will be specified further on. 


TABLE I 
Liquid t у; ç n y, Ay 
33 74.41 0.31 
ы 48 74.63 0.53 
10 10 56 74.16 0.66 
85 75.30 1.20 
Water 

33 69.81 0.31 
hy foe 48 70.05 0.55 
56 70.16 0.66 
85 70.70 1.20 
59 25.48 0.18 
Сус1оһехапе 18 25.80 85 25:18 0.43 
120 26.01 0.71 
Ehtyl alcohol 8 23.82 85 24.16 0.34 
40 30.36 0.16 
10 30.20 56 30.44 0.24 
70 30.52 0.32 

Benzene 
40 27.67 0.16 
30 27.51 56 27.76 0.25 
70 27.83 0.32 
|59 95.34 0,24 
15 25.10 | 85 25.43 0.33 
Batyk alcohol ( 59 24.95 0.24 
20 24.71 \ 85 25.06 0.35 
( 56 33.77 0.14 
20 38.63 l 70 33.83 0.20 
Diosene 56 33.15 0.15 
25 33.00 { 70 33.20 0.20 
33 43.69 0.14 
48 43.73 0.18 
Glycol 47 48.55 59 43.81 0.26 
85 44.00 0.45 
122 44.18 0.63 
Toluene 15 29.09 85 29.60 0.51 
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The value of the r. e. within the ranges a-b and c-d is determined by the composition 
of the liquid and by the rotation frequency. 


Table I brings the values of the surface tension of rotating liquids as compared with 
the rest values, in the ranges a-b and c-d. 
The notation is as follows: 
t, temperature of the liquid in ^C, 
y, rest surface tension in dynes/cm, 
n, rotation frequency in rev/min, 
y, rotational kinetic surface tension in dynes/cm, 
Ay, difference y,-y, in dynes/cm. 
The values for dioxane are presented here for the first time. The values for other liquids 
are recalculated іп dynes/em on the basis of experimental results previously published [2], 
ra. |, [13]. 
The rest values of Table I are in agreement with the values from the literature [12]. 
The values in rotation at different temperatures are higher by about the same amount, 
i. e. parts a-b and c-d probably run parallel to the straight line for rest. 


III. The previous treatment of the liquid 


On the other hand, the width of the ranges a-b and c-e varies with the previous treat- 
ment (thermal and mechanic) of the liquid, as this treatment affects the range of the mini- 
mum b-c, displacing it, or changing its depth. 

In order to investigate the effect of previous thermal treatment, we cooled the liquid 
to different temperatures and subsequently measured the r. e. alternatively at rest and in 
rotation, while heating the liquid. In Table II some values of the temperature are given, 
at which the measurements started, and values at which inversion of the r. e. occurred. 
The values are taken from experiments with glycol [8]. 


TABLE II 
Initial temperature Inversion temperature 
5.45 9.6 
6.35 9.5 
qp V 9.5 
1.10 10.5 
7.90 9.1 
8.25 14.3 
10.90 16.6 
11.15 inversed 
12.20 inversed 


It can be seen that the cooling temperature of the liquid is important for the inversion 
of the r. e. | 
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For investigating the influence of mechanic treatment we proceeded as follows. When 
the r. e. was negative, we shook the liquid energetically (during 3 minutes) and then conti- 
nued the measurements. 


25525 


n= 85 rev/min 


45.13 
10 ft 12 73 14 75 76 77 78 79 20 toc 


Fig. 3. Influence of shaking on the rotational-kinetic surface tension of the glycol: empty circles indicate values 
before shaking, empty triangles after shaking, and S means the moment of the shaking. Full sings represent 
static values of the surface tension in same conditions. 
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Fig. 4. As in fig. 3, for benzene. 


It was found that the r. e. became positive. Following a repeated shaking after about 
half an hour, the т. e. remained positive. Also, within ranges where the r. e. is positive, 


shaking produces no change. 
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In Figs. 3, 4, 5, and 6, the experimental curves are given for glycol!) (Fig. 3), benzene 
[14] (Fig. 4), water [14] (Fig. 5), and butyl aleohol [14] (Fig. 6). Circles denote values pre: 
vious to shaking, and triangle—those subsequent to shaking; S marks the moment of shaking. 


n = 56 rev/min. 


0 
22 15 16 17 18 19 20 24 2 c 
Fig. 5. Аз in fig. 3, for water. 
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ст 


13 14 15 16 7 18 19 £9c 
Fig. 6. As in fig. 3, for butyl alcohol. 
/ 
Subsequent to such mechanical treatment, in order to bring the liquid to its former 
state, it must be kept at rest for several days. Cooling shortens the time required for reco- 
very. 


IV. Other experimental investigations in relation with the minimum 


Another factor influencing the minimum is the time during which the liquid is in rota- 
tion [10]. Fig. 7. shows the surface tension versus temperature curves for glycol, when 
rotational and rest measurements were made alternatively, i. e. when the rotation was 
interrupted regularly. Fig. 8 shows the same curves, though obtained separately, the liquid 
being either at rest or in rotation continually. 


1) These measurements for glycol are new, in agreement with result previously published [10]. 
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It is seen that continuous rotation shifts the minimum towards higher temperatures. 

It should be noted that, for the reproduction of the temperature-dependence curves 
of the r. e., the previous treatment being the same, it is also necessary that the experiments 
be carried out carefully, avoiding sudden heating or cooling and external perturbations, 
such as currents of air in the laboratory. 


n = 85 rev/min 


46.09 ! 
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Fig. 7. Temperature dependence of static (€) and rotational-kinetic (о) surface tension when measurements 
in rest and rotation have been performed alternatively. 


V. Vanishing of the r. e. 


At high temperatures (range е-/), the r.e. was found to be zero. The vanishing point 
depends on the kind of liquid used. Figs. 9, 10, and 11 represent the domain of vanishing 
of the т.е. (range d-e and partly e-f) for different rotation frequencies, in the case of 
water (Fig. 9), benzene (Fig. 10), and dioxane (Fig. 11)?). 

From these curves it can be concluded that the vanishing point of the r.e. is independent 
of the rotation frequency, at least for low frequencies. 


VI. Conclusions 


The foregoing results can be obtained by other methods also. The appearance of the 
r.e. at low temperatures and its disappearance at high temperatures has been demonstrated 
by a bubble method[16]. However, the ring method is the most appropriate, [17], [18], [19] 
and therefore we used it throughout the entire experimental investigation. 

We explained the r.e. as follows [11]. 

Liquids can exist in two different molecular states, one resembling solids and the other 
gases [20], [21]. In the first, the molecules are in some order, constituting relatively large 
groups. By putting the liquid in rotation, a rearrangement is produced, which changes 
the energy state, thus raising the surface tension. This constitutes the positive r.e. 


2) The curves for 56 rev/min in the case of benzene and dioxane, and for 40 rev/min in the case of water 
have been published previously [11], [15]. The other are presented here for the first time. 
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The minimum is due to some special metastable state depending on the temperature; 
rotation frequency, and previous treatment. 

In the second state, the termal agitation prevents any order, and so rotation brings nox 
change and therefore no r.e. appears. 
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THE EFFECT OF LONG RANGE ORDER ON ENERGY BANDS 
IN BINARY ALLOYS 


By А. Совстоукт AND D. CRECU 


Institute for Atomic Physics — Bucharest 
(Received June 20, 1960; revised manuscript received October 23, 1960) 


Perturbation theory is applied for investigating the energy bands of a partially disordered alloy, 
starting from a perfectly ordered one for which the energy bands are assumed to be known, and 
proceeding to consider the perturbation produced by the disorder of the alloy and to compute 
the corrections to the energy levels by perturbation theory in the first and in the second approxima- 

d tion. The width of certain forbidden bands on Brillouin surfaces was found to decrease as the long 
range order parameter, A, diminishes. At complete disorder the respective forbidden bands 
vanish completely. The electron energy is investigated as a function of the propagation vector 
in the neighbourhood of the forbidden band. 


The effect of long range order on the motion of conduction electrons in alloy s is interesting 
not only from the theoretical point of view. Some physical properties of the alloys, especialy 
electrical conductivity, depend essentially on the motion of the conduction electrons in the 
field of the atoms of the partially disordered alloy with definite long range order. The present 
paper deals with the problem of one-electron energy levels of a partially disordered binary 
alloy. We will investigate the dependence of energy bands on the Jong range order parameter R 
of the alloy, and also discuss briefly the effect of the short range order parameter o. 

The problem was first raised by Muto (1938). The case of a one-dimensional alloy 
with various successions of atoms of type 4 and B was dealt with by Landauer and Helland 
(1954), who utilized computing machnies. For completely disordered alloys, interesting 
results were also obtained (Parmenter 1955, 1956; Korringa 1958). The case of a partially 
disordered three dimensional alloy was investigated by Smirnov (1947, 1955), who generalized 
Muto's computations. His treatment concerns especially strongly bound electrons. In 
a previous paper by the present authors (Corciovei and Grecu 1960), preliminary results 
on this problem have been published. 

The great difficulty in treating the problem stated, consists in the absence of special 
periodicity of the partially disordered alloy, notwithstanding the circumstance that the 
atoms of the alloy are assumed to lie at periodic lattice sites. For this reason we will attempt 
to solve the problem by perturbation theory. In the present treatment, the unperturbed 
and periodic potential is characteristically considered to be the potential of a perfectly 
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ordered. alloy, for which the energy bands are assumed to be known (Peierls 1955). One: 
then proceeds to apply the perturbation produced by the disorder of the alloy and to: 
compute the corections to the energy levels by means of perturbation theory, in the first: 
and in the second approximation. 


1. Energy bands in an ordered alloy 


We are dealing with a binary alloy with the components A and В in equal proportion. 
When the order is complete, we suppose that the lattice can be separated into two inter-- 
penetrating sublattices, with an equal number of points, one sublattice occupied only by A | 
atoms and the other occupied only by B atoms. For simplicity we shall deal with a centered | 
cubic alloy. In this case, the two sublattices are: the sublattice of the cube corners 
and that of the cube centers. If the unit vectors parallel to the sides of a cube are denoted | 
by ù, i, d, the elementary cell will be determined by a, = ai, a, = а.а» = ай, where 
a is the side of the cube. 

Let т; denote an arbitrary lattice point of the first sublattice, and r, — an arbitrary 
lattice Palu of the second sublattice. Evidently, we can write 


Tj = mG, + та, + та; T= (mi + 3)G, + (ma + 3) Ay + (mg + Уа; (1) 
where m4, тә, ms are integers. The reciprocal lattice is determined by the vectors 


© „= 2л (n.b, + nb, + n,b;) (2) 


п 


where лу, ng, ng are integers, b,= (a, x a;)lla,(a,>x aj] = (1/а)4,, and similarly b, = 
= (1/а)і,, b,— (1/а) їз. The propagation vectors of the free electron waves are given by 


Әл 
= N (nb, + naba + nsb;) (3) 


where JV? is half of the total number 7 of the alloy atoms. In the free electron approximation 
the wave function of an electron is 


p = exp (ihr) (3) 


which corresponds to the energy E, = (ËJ2m) h?. 
We shall assume that the ОКИ of a single electron in the field of an A ion is U, (r), 
and in the field of an B ion Uy (т). For complete order, the А ions occupy the points Т, 


and the B ions — the points т. Thus, the basic, static, and periodic potential U о(ғ), fol 
the electron in position r, will be 


U, (m) = > U,(r—rj + >) Ug(r— vj) (4) 
j k 


The wave function of an electron in the field of an ordered alloy is obtained by solving 
the Schródinger equation 


(Ho + Uo (r)) vy, = еу, | (5) 
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where H is the Hamiltonian of the free particle. By reason of the periodicity of U(r) 

. . . d 
we can test the wave function as a combination of two free-electron functions (3^) differing 
by a reciprocal vector C, (Peierls 1955), namely 


nis се!" a Cee = суф + соф» (б) 


This yields two solutions for the energy 


Е-Е; E,— E, \? Ч 
a 2 + 5 + [Оса (7) 


the corersponding coefficients being given by 


» » 
cD = | PIE | , = | ned cote : 


=; » HE % 
(оу as al оор Мыш. Vs 1 2 
É | Eid Eg x 4 % | Ib EL x (8) 


where Ug, is the integral 


Uc е ШІ exp (— iG,r) dt (9) 
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О denoting the volume of the elementary cell, and where the phases have been so chosen 
as to obtain the maximum symmetry of the expressions of the coefficients. U,(r) and 
Ug(r) have been assumed to possess the spherical symmetry, whence the value of Ug, 
is assumed real. Moreover, we have chosen the two states фу and @, corresponding to E, 
and E, so that E, >> Е,. Now, in the case of spherical symmetry, we can write for (9) 


D ІШЕ e—:GnR dr + E Ug(R) e- iG» dr 
Q 9 


== U4Gn 5 е—їл(т+п»+Еп) Ова, 
апа 


Do Ї U,(R) sin (RG,) КаК 
ñ 2 


where X = A, B, R is a vector centered at the centre of an A or B atom, respectively, the 
integration being carried out on an elementary cell; G, is the value of the vector ©,, and 
R — the magnitude of the vector R. We shall also assume Ugg, > Ugen: 

The principal result of Peierls’ method consists in the fact that, in the dependence 
of the energy on the propagation vector h, forbidden bands appear, on the so-called Brillouin 
surfaces. A Brillouin surface is determined by the reciprocal vectors ©,, as the normal 
plane through the centre of the respective vector С. In the case of the completely ordered 
alloy, each Brillouin surface, with G, given by (2), is a discontinuity surface. 
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The first Brillouin zone is the cube determined by the planes normal on the smallest. 
reciprocal vectors and passing through the centres of the corresponding vectors. These: 


vectors are 
Eh, doh Eb (10) 


Obviously, the volume of the first Brillouin zone is (27/a)®. 
Let us now consider the following twelve С, vectors 


2л аі), 210), ШЫ) (1) 
The second. Brillouin zone will be situated between the dodecahedron obtained by tracing 
the planes normal on the centres of these vectors and the cube formig the first. Brillouin 
zone its volume also being (2л/а)9. To each Brillouin zone № = /7//2 electronic states are 
known to correspond, if the electronic spin is neglected. If the later is taken into account, 
there аге 4 electronic states in each zone. 

Let n denote the sum 


n= п + Nna + ns (12) 
If there is complete order, the width of the forbidden bands are 


for Brillouin surfaces with odd п. and 


2| UG, + Usa, (139 


for Drillouin surfaces corresponding to ап егеп n. 

In the case of a perfectly disordered alloy we expect that some forbidden bands will 
vanish. In fact, when an alloy is perfectly disordered, any point of the j sublattice is 
equivalent to any point of the k sublattice. In the present case, the vectors of the elementary 
cell will be taken as being (— a, + a, + a3)/2, (a, — a, --а,)/2, (a, + G, — а,)[2, so 
that any point of the whole lattice сап be obtained by adding integer multiples of these 
elementary cell vectors. By writing the reciprocal vectors of this lattice one obtains only 
the vectors G, from (2) corresponding to even n. Thus, the smallest reciprocal vectors are 
those given by (11). 

We can conclude that, in the case of perfect disorder the surface of the first. Brillouin 
zone from the ordered alloys ceases to be a discontinuity surface. The same is true for the 
Brillouin surfaces corresponding to odd n. It is well known that the long range order para- 
meter R varies between 1 (the case of complete order) and 0 (that of complete disorder). 
Thus, we expect that the width of the forbidden bands decreases with R for n odd, vanishing 


at R = 0, and is independent of R for n even. We shall demonstrate this statement in the 
following sections. í 
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2. First order perturbation theory for a partially disordered alloy 
With the aim of taking into account the actual positions of 4 and B atoms in a partially 
disordered alloy, we introduce the parameters ш; and u, defined as follows: 
pai $ A 
Hp if r, contains an atom of type 
|-1 B 
Hg = if r, contains an atom of type | 


(14) 
В 
Thus the potential of an electron in the field of the atoms of the disordered alloy 181) 


U, Gr) = 31 [B (1+ n) Us (r= n) +} (1— n) Us (n— 7)] + 
т > [$ (1— 4) U, (r— т) +} (1+ и)» (r— 79] 


(15) 
The actual potential U,(r) and the basic potential (ғ) from the case of complete 
order differ by the perturbing potential H?*', which takes the following expression: 


нэ = У д) = 0) 
7 
wherein U(r) = (Ú, — 


(16) 

40/2. It should be stressed that, in eq. (16), we have taken into 

account the effect of disordering only and not considered the effect of thermal vibrations 
of the atoms. 


Now, by perturbation theory we can compute the corrections to the energy values 
of eq. (7) in the corresponding states (6), (8). In this section, we shall compute the first 
order corrections. In the next one, we shall discuss second order perturbation theory. 


The energy correction corresponding to the wy state, in the first order approximation 
of perturbation theory, is given by 


1 . 
Ay = 190 I (17) 
INQ 
Clearly, we can write 
ا0‎ 


1 Жы o —iGln" * oiGnr ег 
1740 ШЕГЕ vag | (ee Car 0,026" 587) [Pedy 
2 NO 


z NQ 
and, by eq. (16), the first integral is found to vanish. We now introduce the long range 
order parameter by 


ВЕЕ in 
2 


(18) 


1) Elsewhere, Corciovei and Musa (1960) have considered а more general potential by taking into account 


also the effect of thermal vibrations. One can obtain the potential (15) from the more general potential mentioned 
by equating the thermal displacements to zero. Corciovei and Musa (1960) have improved their potential in the 
problem of the effect of order on resistivity in binary alloys. 
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and, on carrying out the calculations, we have 


1 
AQ = (cica + 163) гуур | Пре exp (— iG) dt 
2 “ 
+70 


ы 3 1 j —in(ni+n,+ ns 
iv (сі + C102) 77 (“ава am Ugan) ps (1 е 0) 6 (mitnstns) У (1 == 2 
7 


== [1 — е- їл(т+т+т)] (Ао, oe Гвс,) 


(cia " C402 )‏ اڪ 


Finally, the correction A™ to the energy levels ғу and e, can be taken in the form 


0122 (Ay = 82) 12 Но £,9)] ^ 1- К ] — e—inm +n +n)| (U 2217 19 
Ау; = 2 FE j eoru 2 [ е ] ( AGn BGn) ( ) 


wherein the expressions of eq. (8) have been taken into account. 
For odd n, expression (19) becomes 


[(Е, e (Е,- £1,9)]^ 
E, + Е, — 286,5 


A®, =2 (1 — R) (Usa, — Unas) (19) 


For even n, the correction (19) vanishes. 


If we consider a Brillouin surface with odd п, the corrections to the energy levels have 
a simple form, namely 


AC — — (1— R) (Ugg, — Unan) (20) 
4@= (К) (U sen т Ugan) 


which shows that the upper level, denoted by e decreases, and the lower level, denoted 
Бу £ rises. The non-corrected values of & and ғә on a Brillouin surface are 


a= B+ Ug.» 8 E— Us, 
wherein 
i Е= Еү= Е, (21) 


and that the corrected values of the energy levels, in the first perturbation approximations, 
are readily proved to be 


s, + 4= E+ ЕШ, 
&+ A — E— RU, (22) 


It is seen that, for complete disorder, the forbidden band for n odd vanishes. On the contrary, 
for even n, the forbidden band remains unchanged. 
From (22) the band separation is 


2RUg, = 2К(О до, Usa) (23) 
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for odd л. For even л, the result is 
206, Te 2(U де is Овал) (23) 


which confirms our previous statements. 
Utilizing eq. (19’) and Peierls’ formula (7), we can write down the corrected energy 
level in the vicinity of a Brillouin surface, in the first order of perturbation. A results, 


for two energies corresponding to two states p, and yy, such that h — h' = 6. (odd n) 
are 
ШӘ G) ы Жж 2 | ТА | Ue 
+41 = E, o us 10” F Uy — | (рген 
Мел e X [(y2-- 091% Р 0% 
2 2 2° + Mos xcd (2+ US (24) 
with the notations 
Е-Е; 


The position of the new level as corrected by the effect of disordering can be discussed 
with respect to that of the energy level in the approximation of free electrons. It is to be 
noted that for 


R = А (26) 


U (0 
the two levels coincide. Conversely, for a given R, no у exists for which the two levels coincide, 
except for R < 0.5. We also note that К, as given by eq. (26), is essentially dependent on 
the difference between the two levels E, and E. | 
In the case of a perfectly disordered alloy, the values become 
edP s E ydg Ure [у-уу] 
в + AD = E,+ у(у*+ Ut) [реч U95— у] (247) 


It can be seen that energy of a level within the zone is greater than in the case of free electrons, 
nevertheless remaining smaller than the free energy corresponding to the Brillouin surface. 
The contrary can be said for the levels outside the Brillouin surface. It is also interesting 
to note that the graph of energy versus the magnitude of the propagation vector normal 
to the respective Brillouin surface has a vanishing slope on the surface. 


3. Second order perturbation theory for а partially disordered alloy 


In principle, the problem of commputing the second order approximation of the energy 
levels of one electron in the field of a partially disordered alloy consists in finding, by 
means of perturbation theory, in second order, the corrections to the energy levels given 
in (24), where the perturbing potential has the expression (16). Unfortunately, the problem 
is highly involved, and we shall be able to solve it only by making some drastic approxima- 


tions, which will allow us to discuss the physical consequences. 
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From the beginning we shall тезїтїсї ourselves to the computation of the second order 
correction of an energy level within the first Brillouin zone. There will be no difficulty then 
in coming to some conclusions about the corrections to the energy levels situated in other: 
Brillouin zones. We may say that an interesting situation arises in the neighbourhood of the: 
first Brillouin surface, as it is the first odd one. 

Let С, denote the vectors (10) determining the first Brillouin surfaces in the case of 
a completely ordered alloy. We shall compute the second order correction of an energy: 
level corresponding to a state y, situated within the first Brillouin zone, with h situated in. 
the vicinity of a Brillouin surface determined by a given G4. In the second order perturbation, . 
we have to consider the matrix elements of the perturbing potential HP* taken between | 
the state y, and the other states yẹ, with h’ > h. We have to consider only, the states yy. 
for which the denominator in the energy correction is small. For this reason, summation. 
will be. carried out only over the states yẹ, with h' situated in the first and in the second. 
Brillouin zones. To avoid misunderstandings, we shall divide this h' subspace into four 
regions. The first region, I, is determined by all the vectors R’ whose projections on G4 have 
the same orientation as the vector G. The second region, II, is determined by all the vectors 
һ-б,, where h' belongs to the first region. Lastly, the third, region, ІП, is given by all 
the vectors —h’, and the fourth, IV, by all the vectors —(h' — G,), where R’ belongs to the 
first region. We shall denote the energy of the free electron, with propagation vector h’, by Ej, 
and that of the free electron, with propagation vector һ'— Су, by Ез. In subsequent 
discussions we shall consider h’ to be a vector of the region I. 

Now, the Peierls energies and wave functions for states situated in the four regions 
can be written as follows: 


Py = сі ей" ТИ с gio Gr 
ї 1) ТАЛ 1 
_ Etk, Dy S uka |. 
PN ao. AIL lc >. + |06, (271) 
кале св" - €; gi - Gor 
f Ey Esq; d d pe aet а 
"h'—G, = 2 % E j ба, (27ш) 
ji od E atti Nro ЕС KN erar 
Yn = Cy € F C9 € 
ЕЕ [Rey ЖЫЛ 
Sec | 5 31061 | -єм (27111) 
EL El та 4 
5 s 1 ш SE "R 9 
S I. = = aa mE 2 АА 
(Nee) grt 2 “| дай | —5»-a, (27 ту) 


where 


22 | e ые ew р ШТ. Қызды bas % 2 
ty + Ey — 26р |” P Ei + E; —2ey | ` f ) 
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Obviously, that the following properties are satisfied: 


* * 
Ven IAS CP tay = w G, 


For calculating the second order correction to the energy level of the state py, we shall 
apply the formula 


Н?Р“, 
AQ ЕЭ pi | | hh 


(29) 


En — Ep’ 


2 per р рег € 
Е |Н Weal | H E П На е, 1 
En — &—һ | 


En — ёһ 6, Eh — 6-һ-о, 


hh 
summation in h' being carried out over the first region I. In view of the summation, we shall 
pass to integration over the first region by means of the formula 


VOR | 
2; ... . е dš h 


Detalis of this integration will appear in the calculation. We now proceed to find the matrix 
elements of HP®. 

We have to compute the matrix elements between the state yy and the states Yp Vy —G, 
Y_n> Y-wie, We note that the state y, and the corresponding Peierls energy is expressed 
in the same way as in (27,), by substituting h’ with h. In each case we obtain a different 
result; all, however, can be written in the general form 


Hy RC ... = Му [2 (1-5 ш) Sp hes > (1— uu) ё] p 
3 Z: 


+ M, р> (1- шека > (1— ша) ei] 


(30) 


(31) 


with K =h’—h in cases I, II, and K =h’+h in cases ІП, IV. The expressions for М, 


and M, in the four cases are 


1 3*7 , 1 * ГА » / 
M, = 177 (cici + c$ сз) Ок; М = V (ci e Ug a, + cà 1 Ux a) (321) 
1 [7j 5 7 WE: 1 ги * 7 
Mi = 177) (с1 Ap КС? c2) Uk; M, = 177 (с1 Сә, Окс, == Шум! Ок+в) (3211) 
1 ry * 4 1 * * 
M, 5E 1171) (с\ Cy Он Я7 Сэ, Сә UxK-24,); M, = 400 (сі cà - Ca ci) к-с, (32) 
1 * * ci * Jf + 
M, = 1⁄4 (ісі Ок — с с Ug-2a.); М = 377 (сі c2 — cà ci) Un-a, (321v) 
wherein we have introduced the notations 
Ок-ю, = E | U(R)dr, £ 2-5 0, + 15 2 (33) 


with K = h/ — h in cases I, IL, and K =h’ +h in cases Ш, IV. We note that, in the hypo- 
thesis of spherical symmetry of U(r) the expressions (33) are real. This is not a significant 
limitation, and the transcription of the various formulae in the complex case is evident. 


X 
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We сап now compute the square of the modulus of the matrix elements. The general 


result is 


IH]? = (M, + М) 3) (1— n) — pj) r + 


JJ 


+ (M,— №)? У (1 — n 1— a es n) — (34) 


2 
> 


— (М1— M3) Xi 1.— p) ule + етеу") 
58 
where М, and M, are those of eq. (32), and the same convention is made concerning K as 
in eq. (31). 

Summation in the parameters и, and jj is not simple. In another paper, one of the 
present authors (Corciovei 1959) has discussed this problem in more detail in connection 
with the effect of order on electrical resistivity. In fact, the formulae of that paper differ 
from the ones appearing here; however, to circumvent the difficulties arising, we can carry 
out the same averaging as in the earlier paper. 

Firstly we remark that, from (14), и? =1, и? = 1. Thus, we can write (1— uj) = 
= 2 (1— uj), and similarly for k. The long range order parameter R has been introduced 
in (18). The short range order parameter с, for two neighbours j and £ of the first гап, 
can be introduced by the formula 


o=— ) шш (35) 


where Q is the number of neighbours of first rank in the alloy, and the summation (J, k) 

is carried out over all pairs of neighbours of first rank. We note that it is also possible to intro- 

duce other short range order parameters for neighbours of an order exceeding unity (Gibson 

1956). In the present paper we shall consider only the parameter с introduced in (35). 
Averaging in eqs. (18) and (35) we can write 


~ = — 
ш =U, = Е ци =o (36) 
` چ‎ . . 
where, іп Hj,» J and k are first rank neighbours. We assume the parameters Mj, р, М, 
Шу for neighbours of order exceeding unity to be independent, whence 
=т= зы» ур ص‎ | кору мл Л 
Буй“ Шуну = Ht = Hle = Шш, = (уй, = R? (37) 


where j 5 j’, k z^ k' and j and k are neighbours of a rank > 1. Now, it is readily proved 
that 
-p = > m= (1— RD + (1— R)? 


— 


“(= (0.— =(e— R44 (1— А)? (38) 


2) The neighbours of first rank are the nearest neigbours in the alloy. 
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where j and £ are neighbours of the first rank, and also that 


usns =O), ODO 0-4) 


(1— 4) (L— д) = (1— R)? (38? 
where j and k are neighbours of а rank > 1, and ;Z J, k zz k. 


With the foregoing formulae, we can obtain the average value of (36), Dropping the 
symbol of the average, we have 


Hal =A [(M 1+ M2) (L— R?) — zfx (М1— M?) (c — RD] + 


+ (1 — R)? КМ, + M,)? 2 EEr) | (M,— M)? 22 eir) 


кк” 


ее (М1- M3) p [E 4 gather dp (39) 


where the convention of eq. (31) has been made concerning K, and where we have introd- 
uced the notation 


fk = E Y e+iED (40) 


z representing the number of nearest neighbours of an atom, and summation being carried 
out over all the nearest neighbours of an atom, D being the vector connecting the atom 
with its nearest neighbours. 

Before proceeding further, we shall make the assumption that the quantities U, are, 
independent of К and constant. This assumption, called the isotropic aproximation, is 
currently applied in the problem of resistivity of alloys (Gibson 1956). In the theory of 
resistivity, the isotropic approximation provides a way of finding the significant features 
of the phenomena. We shall make this assumption here in view of avoiding an excessive 
amount of computation. The physical characteristics of the phenomena will not be obscured 
by the isotropic approximation. Concerning the expression (39), it should be said that the sum- 
mations > exp (К?) and > exp (iKr,) vanish, and thus the expression in the great brack- 


ets in eq. * (39) vanishes. А i case exists, namely that of K = 0, when these summa- 
tions do not vanish. This case can occur only for the region II and corresponds to the condi- 
tion h' =h 3). In this single case, the expression in great brackets in eq. (39) does not 
vanish, but provides his contribution in (29), which must not be integrated in eq. (30). 
Let us consider the normality conditions for the coefficients of the wave functions, which 


yield 


, , , , 72 
2 | Гэв 
(сус, E Coco) + (сүсә, + сос)? = 1 + 4сүса 616, 


(сусу = 06) — (суса + C261)? = (с — ch) (сү — 4) (41) 


3) The condition h” = h cannot occur in region I, for the nondiagonal matrix elements. In regions HI 
and IV, K = h’ + h and therefore does not vanish. 
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In this manner, we obtain the following expression for the averages of the squares of ` 


the matrix elements: 


20% 
30. 


ЇЕЭР" = 


emcee ST с. 


(1- R2) (1--4сүу4) — zfx(ct — с) (^ — ca) (o — R?)) — (421) 


20% 
100 


|А? = 


[(1 — А2) (1 — 4с,сосусо) — zf к(с — о) (ci? — es?) (с — RI] + 


+ 4U&(1— R)? T rp (42) 
2 
н wl? = FUE qa je (1+ asd) — ау} — d) (43-40) (o — RD) 0249 
er ын Ё Т 
НЫ —w+e,|? = К?) (1 — 4сусьсусь) — zfk(c$ — c3) (ei? — cs?) (o — R?)] 


(421v) 
where K = h' — h in (42у) and (425), and K = h' +h in (42,1) and (42,,); y and U 
are given by (25). 
It is now possible to write down the complete expression for the second order correc- 
tion of an energy level corresponding to a state y, in the first Brillouin zone, region I 


2 
4j = PO — RD + Qi — Жу -2Uk0 R 49) 


where 


1 
— t 
E, — EV? % ГЕ — ЕЗ % (1) , 
( = 1 ШЫ | 5 | + [Walt] SN EN (43) 


Uk Ep re 
| DRRR, E 5 x TUS Л : (487) 
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the notations 


Їл Буд). уаш 2 11% 
ашы ані; ü ( ^g | ШИ 


, 1 7 ЯМ чүл? 2 
в? = БІЗ к 4 (8 5 zi + |Ug, 


2 

5 i Ш 
being also used. It should Бе noted that the last term іп ед. (43) arises from the great brac- 
kets in (39) and that in formulae (43’) and (437), which derive from (29) and (30), integra- 
tion is in region I and we must consider that the poles are eliminated. Thus, in eq. (43), 
the significant contribution is obtained by integrating in the neighbourhood of h' = h. 
In order to carry out the integration in the neighbourhood of poles, we have to take the 
principal value of the integral. 

We shall now compute the term given in eq. (43) only. This term provides a signifi- 
cant contribution in (43). On the contrary, the term from eq. (437) contributes insigni- 
ficantly to (43), especially in the neighbourhood of the Brillouin surface. Indeed, in (43) 
this term is multiplied by (ø — А?), a quantity which is negligible in comparison with 
(1 — R?) or (1 — R)?, especially for small R. On the other hand, eq. (43”) contains the 
factor (E, — Е) which also tends to reduce the contribution of this term as compared with 
the other terms іп eq. (43). We note that factor (E, — E,) allows the vanishing of (43") 
on the Brillouin surface. Lastly, eq. (43") contains the factor fg, which also diminishes 
the value of this term. 

Thus, one is justified in dropping the term containing the factor (с — R?) in eq. (43). 
We shall now consider the first term іп A, namely (1 — R?) P. As we have already pointed 
out, the significant contribution in this term arises from the principal value. Conversely, 
when h’ is not in the vicinity of В, the integrand in Р is very small. It is a usual approxima- 
tion (Parmenter 1955) to assume the contribution on the surface h’G, = 0 of the integral 
(437) to vanish. When № is in the vicinity of the Brillouin surface, the contribution for h’ 
on the surface is significant and will be taken into consideration. 

We shall compute P for h situated in the neighbourhood of the Brillouin surface. 
In this neighbourhood we can introduce the Peierls effective mass, in view of avoiding 
excessive computation. This is also a common approximation in such theories (Paramenter 
1955). The introduction of effective mass is equivalent to an expansion in series in Peierls’ 
formula for the energy. Namely, in the neighbourhood of the Brillouin surface, we can 


write for (7) 


h? G? Mi? h*/2m x9 2 ЗУН 
(LA с=т= у= — —— ‘os == БЕ 
Un 2m 4 + |С E l + 5102) 22 : а 
n* 1 2 р 
+ — (1+ — ЖС соз? O | А"? (7) 
2m 4 


where Ө is the angle between h” = (G,/2) — h' and G}. This approximation will also be ap- 
plied for eq. (44). 
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In order to carry out the integration, we shall shift the origin of coordinates to the 
point determined by 6,/2 and make use of spherical coordinates in h space. We have to 
be careful to take the principal value correctly. It is not necessary to insist on this point 
here. The calculation is not difficult, but very lengthy.*) The final result is: 


G, -n | 
PE ee ; 49) 
(21) * Es: W (y2 + 02)% Ев. IV |G, : 
Ны ежен T 


We note that in this expression, the essential contribution arises from the first term inr 
eq. (437). pw 

Finally, the energy level, corrected in the first and in the second approximation, fort 
a state Yp with h situated in the first Brillouin zone, becomes 


I UA U? 
--42-42-28, + [упр Upa ЛЖ U4 — у]— R sr pays t AN (46), 
where 
Q а U 1 
Саре “үлі. : 1— К?) — 
Ар 242 Uk Esa W T 0941 Ea WG, _ - ( ) 
U 2 
2U$ dete N EER (47) 


(у + U3) 


What are the principal conclusions from the preceding calculation? In fact, the funda: 
mental result is that the second approximation of the perturbation theory does not signifi- 
cantly affect the first approximation, which remains the predominant one. Indeed, it is an. 
important result that 49) vanishes on the Brillouin surface, and that its dependence on the: 
magnitude of the propagation vector normal to the respective surface presents a vanishing: 
slope there. Thus, on the first Brillouin surface, the effect of the first order perturbation 
remains unchanged, namely the disappearance of the forbidden band, in the case of complete: 
disorder, and also the vanishing slope of the energy dependence on 2n - l. At the: 

1 
same time the band width for the first Brillouin surface remains, even in the second approxi- 
mation, namely 


2R(U sa, — Ово.) (49) 


In general we expect that the second approximation does not modify the results (23) and 
(237), which give the band widths, in the case of odd and even n, respectively. It should be 


^) We are interested in the significant feutares of the result. Thus we have applied, instéad of (7^), the 
simple formula ey, = Ep(1 + W cos? Ө 1/2) which represents a good approximation of (77). 
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stressed that the foregoing results are still valid even if we were to take into consideration 
the term (o — R?) - Q in eq. (43). 

Fig. 1 shows the one-electron energy dependence on the propagation vector h normal 
to the first Brillouin surface, in the neighbourhood of the surface. We have plotted the free 


Brillouin 


energy 


energy of free 
electrons 


energy for complete 
order 


energy in first 
approximation 


energy in second 
approximation 


propagation 


vector h 


Fig. 1 


electron energy, the energy of a completely ordered alloy, and the corrected energy in the 
first and іп the second approximation, for the case of a perfectly disordered alloy. The 
second approximation is seen to yield a negative correction within the Brillouin surface 


and a positive one outside. 
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STATISTICAL THEORY OF MULTIPLE PARTICLE 
PRODUCTION WITH ANGULAR MOMENTUM CONSERVATION* 


By Zino Kopa** 
Institute of Nuclear Research, Polish Academy of Sciences, Warsaw, Poland 
(Received September 17, 1960) 


A new formulation of the statistical theory of multiple particle production is investigated 
in which the angular momentum of the system is strictly conserved, and which can thus be 
applied to the analysis of peripheral collisions and the angular distribution of secondary particles. 
Fundamental formulae (§ 2), simple examples (§ 3), and some relevant remarks (§ 4) are presented. 


§ 1. Introduction 


Since the pioneer work of Fermi (1950) the statistical theory of multiple particle pro- 
duction has been repeatedly discussed and applied to a number of problems, and various 
modifications, improvements, and refinements have been proposed!. 

As far as I know, however, no serious attempt has so far been made in order to take 
into account the angular momentum conservation, except the thermodynamical treatment 
of Fermi (1951) himself?. 

At present we know on the other hand that there are many cases in which large angular 
momentum transfer in the collision plays an essential role?. 

In the energy region of several GeV, for example, a considerable fraction of the events 
are regarded as due to peripheral collisions. The sharp anisotropy of the angular distribu- 
tion of secondary particles in the center of mass system of cosmic ray jet showers indicates 
that at superhigh energy a large angular momentum must be carried away by the assembly 


of created particles. 


*) A brief summary of this article has appeared in Nuovo Cimento 18 608 (1960). 

**) On leave from the Research Institute for Fundamental Physics, Kyoto University, Japan. 

1) There are quite a large number of works devoted to this theory and its variants. The readers may be 
referred, for example, to some review articles, where more extensive lists of literature can be found. (Rozental 
and Chernavskii 1954; Milburn 1955; Belen'kii et al. 1957; Ito and Tanaka 1958; Koba and Takagi 1959; 
Hagedorn 1960). 

2) This approximation can be valid only for a large number of particles at superhigh energy. For more 
details concerning this point, see $$ 4.3-4.5. 

3) See § 4.14.2. 
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In this note I should like to discuss a new formulation of Fermi's statistical treatment | 
of multiple particle production, in which the condition of the angular momentum conser- - 
vation is explicitly imposed and which can thus be applied to the study of peripheral colli- - 
sions and the analysis of angular distribution of secondary particles. It is hoped that such 
an approach will be useful for a more unified treatment of head-on and long-distance colli- - 


sions. 


$ 2. Formulation* 


Let us discuss the simplest case where all the secondary particles: have the same mass : 
and zero spin, and no other degrees of freedom (neutral pions, for example)?, since the : 
generalization to more complicated cases is straightforward. The conventional formula for : 
evaluating the relative magnitude of the phase-space volume available for № particles with i 
total energy W and total linear momentum Р 156 


_ dUy(W,P) 
where U is given by 
: ПО ЗЭН а тр 
= SIG i Зу. 
шэг A fon) 
ка], со 
x ó(P — Ур) Ө(Р- У ш), (2.2)! 
with | 
[gx E R3: effective interaction volume’, (2.3) | 
w; = Vp? + и?: energy of secondary particle. (2.4) | 
O is the step function defined by 
— ПИ O 


The factor 2 M/E, M being the mass of the incident nucleon and E the total energy in the 
center of mass system, represents the effects of the relativistic reduction of the interac. 


3) The natural unit system 4; = c = 1 is used throughout this note. 

>) Experimental data seem to suggest that we should treat the nucleons on a different footing. So we are 
concerned here with the equilibrium of pions only, assuming definite values of energy and momentum — 
later, angular momentum instead of linear momentum — transferred to the pions. 

6) Usually P is put equal to zero. 

7) Usually R is taken as l/u, the pion Compton wave length. 

8) We do not identify E with W + 90), Thus ТЕ — 2M) (inelasticity) is an adjustable. parameter. 
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tion, which is usually visualized as the Lorentz contraction of the effective volume?. 
The factor 1/N! is due to the indistinguishability of secondary bosons. 
First we notice that (2.2) can be transformed into 


d. тығады ° 
Un(W, P) = CaN NT TT ШЕГІ IV; 21 X 
pail e —oo 


N 
х |Пехр (рад 0P — X p) 6 (W — У) w). (2.6) 


1-1 
The Lorentz factor 29R/E is replaced here by a more general function €; which can also 


depend on ae; and p;. A more covariant version of the statistical treatment is obtained, for 
example, by putting!? 


g = 4 (2.7) 


It is possible, further, to introduce an arbitrary form-factor (non-uniform interaction vol- 
ume). See § 4.2. It is to be emphasized that the expression (2.6) describes explicitly the 
original picture of Fermi (Fermi 1950; Ishida 1959): the square of the transition matrix 
element is regarded as being proportional to the probability of finding simultaneously all 
the secondary particles, each of which is represented by a plane wave normalized in a large 
unit volume, in the small effective volume of interaction. 

In (2.6) the final state of the system is specified in the representation in which each 
particle is labelled by its linear momentum p, and the total linear momentum Р is a good 
quantum number. Now we can use, however, an alternative representation, in which each 
particle is labelled by its angular momentum l and one of its components m together with 
the magnitude of linear momentum p, and further go over to a representation in which the 
total angular momentum of the system is a good quantum number. h 

For this purpose we remove the factor 0(P— Хр) from expression (2.6), expand the 
plane wave into spherical one, 


со 1 

exp (ipx) = т> Pr ДОШ) Т, (а), (2.8) 
=0 = 

where w and o, represent the solid angle in the coordinate and momentum space respec- 
tively, integrate then with respect to œ, and obtain, the assembly of the product of one- 
particle eigenfunctions of angular momentum. This corresponds to a representation where 
all l's and m’s are diagonal. Wé make a further unitary transformation by taking suitable 
linear combination of them, so that the total angular momentum becomes diagonal. ‘Thus 


we construct eigenstates of total angular momentum and restrict ourselves to a particular 


э) We cannot take this “contraction” literally. If the effective interaction volume is really contracted in 
the center of mass system where two incident particles stop (in the case of inelasticity one), the angular momentum 
conservation would be violated. It should be noticed that in the usual treatment the so-called “contraction” 


simply means multiplication by a reduction factor 29X/F. 
10) Such a version was first discussed by H. Fukuda (1953), and later by Ito and Tanaka (1958); more 


recently Srivastava and Sudarshan (1958) developed a practicable treatment starting from this assumption and 


making full use of relativistic covariance. 
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one of the latter. In this way the relative value of the space volume available for V particles with 
total energy W and total angular momentum L, one of its components being M. will be 


given by the following expression! : 


dUnx(W; L, M 
Syl IV; ГА М) = — 4 E , (2.9) 


(L) 1 


Li gea IET cte 
А БҮЙ Пу! 


28 ч 
Ux (IF; L, M) = [ š) 
IE т РЫА 


м R Ex 
X FH [f r; dr, Гао, | p;dp;S, (Е, W, Pis 21 X 
i=1 0 0 


х [Уу Ky (ШМ; L; Ay, «bya Anas In м| mi іт) Х 


Phys HIN 
X Yim, (01) -+ Yinmy (ON) J, (Рат) --- м (pyra) Be(r— "y ШАР (2.10) 


where лу, nə, ... are the numbers of particles with identical angular momentum. In other 


words, 


and 


Lati aaa 7l, etc. 


To TCU ‘ 5 ati y (D) дүүс 
In the expression (2.10) the first summation U Oro las <s IN: А... aw Q means that all 
possible combinations of integer values for L, la, ..., ly. Ay. Aas Ауа which fulfil the 


following triangle relations should be taken: 


l — 4| SL <S + 2, 


„01.2.6039. 11208 


My BESO ләә ТЭХ ТУУЛ. ЭЭ а иене ales ЭЭГ? 


For. a xed gelo ЫТЫ ЫРА Б Амо: the second summation goes over all 
values of (STS 110,00, М). 


Ihe coefficients Ky(L, M; 1, À,, ... Ама, х» ly] lh, mis іуту) are the matrix ele- 
ments of the unitary transformation from the representation in which lis тү, ... ly, and my 


1N Н ? e š : : 
) Неге we use the wave functions for free particles and thus neglect secondary interactions among the 


produced particles, just in the same way as in Fermi's original treatment. The effect of secondarv interactions 
must be studied separately. A 
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are diagonal, to the one, in which J, ... ly, А, ... Ау_ аз well ав L and M are diagonal. 
Тһеу can be easily constructed by the formula: 

R CEM a las RT IS Л ha АЕ Б.д so Ly, rss) 


-= > КОШ Dy АА та» А. qa) Ky; bon Ags coy 


кл у с та, озо у ys) 
= >) K,(L, M; L, lh, m; A, My) Ko(Aq, мл; 15, А la, тэ: А, ug) X 


E AN 5 
X K,(Às, шь; ls, А3113, тз; Ag, Мз)... 
КЫА, Им-а3 Incas Islas ™n—13 I my), (2.12) 
апа 
IAU, l ты: а, то) CL. ту mas М); (2.13) 


where С (L, la, L; mı, тә, M) is the usual Clebsch-Gordan coefficient. From the well- 
known properties of the Clebsch-Gordan coefficients it is seen that in the multiple summa- 
tion in the right-hand side of (2.12) only those terms remain for which u, = M — my, 


М-3 
ә = M — m, — mg, ..., uy, = M — >} т. Consequently, 
3-1 
KL, My ty, Ay, AN a E IN |» та:2--: Ly, my) =O, 
if m тә. O М. (2.14) 
Also we have from the unitarity condition 
> ТКМ. М; А, А, ... Аз, іма I| hs та; ---3 De my) |? = 1. (2.15) 


mim +... Tm N= М 
We can perform the integration with respect to the angles œ; in (2.10) and with the 
help of the orthonormal relations 


| deo YZ, (o) Yu, (9) = б, (2.16) 


and (2.15), the result is simplified into 


N 1 
Unx(W; І, М) = E ТЕ peri a pb. 


л dme nb Лу! 
N R oo 
Ж Г] [f dr, | ар; (E, W, ри) { рт (pr) OW — ОХ), (2.17) 
22 0 0 ps 
Bri: Ë 
2\N (L) 1 
Su(W; 1, M) = Б DE AL, hs- Iw) TT 
L... IN 
N R со 
х ПИ dr, | dp,&; (E, W, pr) (bz; (pi) 9 (W = 51%). (2.18) 
= © 9 Е 


12) Diferentiation of €; with respect to W is here neglected. 
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Here we have replaced the summation over the parameters 21, ... Ay. ә by the coefficient t 
Z (L; l, ..., ly), since the summand does not involve A's. The coefficient Z (L; Hh, ... 1,)) 
represents the multiplicity of the state with total angular momentum L, M, involved in the 4 
combination of the one-particle states with Jy, ls, ..., ly. 

In (2.18) the factor {pr J,( pr)}® is appreciably different from zero only when pr > 1,, 
because it behaves as!? 


х21-2 


{ай(х)}* == 1011) ЇР for x < I. (2.19) | 


It plays therefore the role of cutting-off small values of p when / is given, or that of cutting- - 
off large values of / when p is given. Evidently this factor describes the circumstance that : 
for lower energy the partial waves with higher angular momentum cannot penetrate into the > 
effective interaction region. In contrast to the case of the usual plane-wave representation, , 
our interaction volume is non-uniform, due to the presence of the angular momentum bar- - 
riers. 

Formula (2.18) gives the relative probability of finding N secondary particles when i 
the total energy and the total angular momentum transferred to the assembly of these parti- : 
cles are specified. In this way we obtain the multiplicity distribution. It depends only on: 
the magnitude of L and not on the value of its projection, M, as could be expected. 

The angular distribution of the secondary particles, on the other hand, depends explic- - 
itly on the azimuthal quantum number M. This we can see by leaving one angular inte- - 
gration (with respect to @,, say) and one summation over 4 unperformed when we go from | 
(2.10) to (2.17) and (2.18), and thus by expressing Sy in the following manner: 


SAFEL M) = Гао M YQ) FD (W,L,M), (2.20) | 
lm 
with 
1 92 x40 
ON) ыл ус ^ 2 = ] 2 

Fu 21 IC(L 2, L; m, M — m, М) x 
R œ% 

Х | dr dp (E, W; p, r) Up pr)? Sy 40V — w; A, и). (2.21) 

0 0 


Here the summation is taken over all values of 2 which satisfy the triangle relation together 


with L and /. Тһе average relative (unnormalized) angular distribution А(о) can be found 
by summing over all possible multiplicities : 


А(о) = 21 | Vig (o) SY FO. (2.22) 
lm N 


In deriving expression (2.21), we have made use of the second line of (2.12) as well as of the 


fact that Sy. , (IV— w; 2, и) does not depend on џи. 


19) ji) -/ БЭЛ 
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Š 3. Simple examples 


In order to obtain some results which can be directly compared with experimental 
data, it is necessary to develop a reasonable approixmation method of estimating the inte- 
grals and the sums involved in expression (2.18), just as it has been done in the ordinary 
statistical theory of multiple production. 

In this section, however, let us restrict ourselves to simple examples of relatively small 
available energy, where we make use of straightforward evaluation for the summation and 
a rather crude approximation for the integration. The errors introduced by such a procedure 
have not yet been estimated!4, Consequently the numerical values given below should not 
be taken too seriously. They are rather meant for illustrating some of the characteristic 
features of our approach. 

The examples we examine here are: i) W = би, К = 1/u, 5, = const = л/6, and ii) 
W = 10u, R = l/u, <, = const. = л/8. If we interpret the values of the parameter € as 
really representing the Lorentz contraction of the nucleon-nucleon interaction volume, 
they should be identified with 29/7, so that the above sets of numerical values would 
roughly correspond to the nucleon-nucleon collision of: i) primary kinetic energy (in the 
laboratory system) 5 Gev, inelasticity 0.5, and ii) primary kinetic energy (in the laboratory 
system) 11 Gev, inelasticity 0.5. The multiplicity distribution and the distribution of the 
angular momenta of secondary particles!? will be worked out for all the permissible values 
of the total angular momentum transferred to the assembly of the secondary particles 
(Lmax iS 5 in the first case and 9 in the second). Неге we take the point of view that the 
inelasticity and total angular momentum of the secondary particles are not necessarily 
correlated. See 8 4.3. 

The first simplification we introduce is to replace the “cut-off” factor {prj,(pr)}? in 
(2.18) by a step-function gO (pr — 1), where g is a constant of order one 16 and © is defined 
by (2.5). This approximation corresponds, roughly speaking, to the classical treatment 


of the impact parameter. Then we have 


у На НОНЕ (322 
Lo... IN 


where 


л 


осе 21 1 
sN (АҺ, ls, oade Ix) = ES nione = x 


x Га ЛЕЛ (р; — 1)] exp ia(W — >j w;) 
-со-ів 14-10 
s AEN, Tr 
= i da exp (iaW) I], (3.2) 


п. 27 xm 
—со—1& 


14) See, however, remarks їп a mote by Z. Koba and М. Majewski, which will be published 
shortly. 

15) This gives a rough idea of the angular distribution. 

16) In the following numerical caleulation g is put equal to one. 
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with 


tg а j: [ө — 1) exp (— iaw) 


| 259 m i 
= р (в 4) exp (— їшї). (3.3) 


IIR 


An approximate evaluation procedure for J, is given in Appendix 1. The results are 
as follows. 
For / = 0, i)extremely relativistic approximation, 


2 1 
ар (3.4) 


and ii) non-relativistic approximation, 


2 "m ; 
Тас V5 «екү/ É exp (— iau). (3.5) 


For 1 > 1 (always extremely relativistic), 


l 1 1 
L == т %Кер(- aR] $e 
Rj | ix iat =] (3.6) 
From (3.2) —(3.6) we obtain the relative magnitude of the phase space volume for the 
allowed configurations (Һ, lg, ..., ly). (For more details see Appendix 2). In Table 1 the 
TABLE 1 


Relative magnitude of phase space volumes for the case RI = 6, € = 2/6 
AL; (he 1» .... IN) X SN: [Ae .... ly) 


Configuration L=0 L=1 L= 2 L=3 L=4 L=5 


(0) 100 da n Фт 
(1) — 99 T = 2% 
(2) -- =: 86 2 25 
(3) = кет = 63 
(4) 2m = 
(5) = - — - 18 
(0,0) 100 = Ж 
(0,1) — 134 25 = 
(0,2) яг 5 62 
(0,3) € = = 30 
(0,4) 22 
(1,1) 35 35 35 


Үүл 
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TABLE 1. (Continued) 


(0,2,2) 0 
(1,131) 3 
(1,1,2) 0 
(1,1,3) = 
(1,2,2) 0 

‚ (0,0,0,0) 1 = = 2 E. Е 
(0,0,0,1) = 1 - e ГА M 
(0,0,0,2) = 7» 
(0,0,1,1) 0 0 
(0,0,1,2) s 0 
0 

0 


(0,1,1,1) 0 
(0,1,1,2) 0 
(1,1,1,) 0 0 
(1,1,1,2) 0 
(0,0,0,0,0) 0 
(0,0,0,0,1) = 0 с - 23 LES 


(0,0,0,1,1) 0 0 0 зэ за S 


¿ 77 means that the contribution vanishes exactly, due to the zero value of Z (L; D, while “0” means that 
the contribution is finite but is negligible within the accuracy of our numerical estimation. 

'Those configurations which are forbidden by more exact expressions for the energy (w — yp + и?) 
are excluded. (For example the configuration (00111) is allowed in the extremely relativistic approximation, 


but it is not listed here) 
Normalized as s, (0) — 100. 


numerical values of these phase space volumes sy(l,..../y) multiplied by the weight 
factor!” Z(L; l,, ..., Ly) for case i) (W = би) are listed. With the help of this table we can 


17) This weight factor gives the number of states with total angular momentum Z and one of its components 
M involved in the given configuration (l, ls, ..., ім). In other words the degeneracy factor (2L + 1) is not 


included. 
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easily find the multiplicity distribution and the distribution of angular momenta of the 
individual secondary particles. The results are summarized in Tables 2 and 3. 

The corresponding distributions for case ii) (W = 10u) are given in Tables 4 and 5. 
The list of phase space volumes multiplied by the weight factor Z(L; h. ..., Ly) is rather 
extensive and is shown іп Appendix 318, 

In working out these values the extremely relativistic approximation has been used 
for particles with 7 > 1. For the S-wave particles we have applied the extremely relativistic 
approximation in some cases and the non-relativistic approximation in other cases, depending 


TABLE 2. 
Multiplicity distribution of secondary pions for the case RIV = 6, Q = л/6 


N 1510 185511 L = 2 Ше) | L=4 | L = S 
1 38% 29% 40% 43% 64% 92% 
2 52% 58% 41% 49% 34% 7% 
3 9% 12% Юс 8% до 07198 
4 0.5% 0.6% 0.4% 0:195 0.095 0.095 
Б 0.0% 0.0% 0.0% эн 2s x 
N 1.72 1.85 1.75 1 6 | 138 ШЕТ 
ТАВГЕ 3. 
Distribution of angular momenta of secondary pions for the case RW = 6, Q9 = л/б 

[ 125540) ТА HolL m L=4 | L=5 
0 77% 30% 24% 16% % = 
1 2295 64% 40% 23% 14% 6% 

2 1% 6% 51% 17% 7% 2% 

3 ы 0.0% 396 4496 1396 376 
4 A Эн E 0.495 S0% 295 

5 — — -- - - 85% 

1 0.24 0.77 | 151 | 1.90 | 2.93 | 4.60 

TABLE 4. 


Multiplicity distribution of secondary pions for the case КІР = 10, 9-018 


N пОИалгагсагтагагтагсртгт 
1 | 18% | 10% | 10% | 11% | 14% | 19% | 29% | 47% 67% | 100%, 
2 | 44% | 46% | 48% | 51% | 54% | 57% | 58% | 49% 33% 0% 
3 29% | 34% | 34% | 32% | 28% | 22% | 14% 4% 0% 0% 
4 7% 9% 8% 6% 4% 2% 0% 0% 0% 0% 
295 096 0% 0% 0% 0% 0% 0% 0% 0% 0% 
N | 2.3 | 2.4 | 24 2.3 23 2.1 | 1.8 | 1.6 | 1.3 | 1.0 


15) This corresponds to Table 1 of сазе i). 
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ТАВГЕ 5. 
Distribution of angular momenta of secondary pions for the case RW = 109 — л/8 

1 |1-0|1-1| і-2|1-з|1-4|1-5|1-6|1-7|1-8| 1-9 
0 5570 30% 2595 2195 18% 16% 159 1105 7% — 

il 5200 48% 84%, 28% 2295 18% 14% 8% UO 0% 
2 ПОА 179% 31% 2999 19% 16% 129% 9% 596 0% 
3 206 4% 8% 299 16% 14% 11% T% 5% 0% 
4 0 0% 1% 595 2195 1395 10% 896 595 096 
5 — 095 0% 1% 395 23% 11% 995 596 0% 
6 — -- -- QU 096 1% 2005 UG 5% 0% 
7 -- -- -- — — 096 1% 42% T% 0% 
8 — — — — — — — 0% 57% 0% 
9 = == = — — — — — — 10095 
1 | 0.6 | 0.9 | 119 | LZ | 2.1 | 2.1 | 3.4 | 4.6 | 6.1 | 9.0 


on the amount of kinetic energy available. The former cases are denoted by 0 (without bar) 
and the latter by O (with a bar). In some intermediary cases we have taken the average of 


these two approximations; such cases are denoted by 0 (with a circumflex accent). 


$ 4. Discussion 


1. We hope that the approach proposed in this note can be applied to various kinds 
of high energy nuclear interactions. It seems more reasonable, for instance, to classify , 
the states of the final pions emerging from the low energy nucleon-antinucleon annihilation 
in terms of angular momentum eigenstates, because then we can make full use of the selection 
rules. 

One of the most sensible application will be concerned with so-called peripheral nucleon- 
nucleon or pion-nucleon collisions at accelerator and cosmic-ray energy regions. Of course 
it would be too optimistic to imagine that such a simple statistical consideration could 
suffice to explain the detailed aspects of high energy nuclear collision. Possibly a more 
specific mechanism, such as the formation of an excited nucleon or meson-cloud (isobars 
or “fireballs’’) will have to be introduced. We shall be more justified in doing so, however, 
when we have eliminated the effects of the “kinematical”? factor due to the relative magni- 
tude of available phase space volumes. 

2. It may be remarked that the concept of a “peripheral? or “long-distance” collision 
is not quite clear-cut. A reasonable quantum-mechanical definition for this kind of collision 
is given by specifying a relatively large value of angular momentum, in comparison with 
energy transferred from the incident particles to the assembly of secondary particles. Obviously 
such a case can be fully studied with our formulation. 

An alternative way, which is much more classical, of defining peripheral collisions 
will be to introduce a particular interaction volume which is effective only at the boundary 
region (interaction volume of a spherical shell type). Such a treatment can be also realized 
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within the framework of our formalism by modifying the reduction factor € into a function 
of the radial distance r. (See the remark after (2.7)). 

Further we notice that results of our formulation are not identical with those of the 
conventional statistical theory, even when we are interested in the average behaviour of 
nuclear interaction, including peripheral as well as head-on collisions. Then, we have to 
take, for example, the weighted mean of our results for specified 1/5 with the weight (2L -- 1) 
up to the maximum allowed value of L. The average multiplicity distribution thus obtained 
will be more or less similar to the prediction of the conventional statistical theory in terms 
of plane-wave eigenstates. But our formalism can afford a further information on the average 
angular distribution of the secondary particles, while this is not the case with the usual 
treatment, as has been emphasized by Hagedorn (1960). 

3. We have mentioned already that the problem studied in this note was once discussed 
by Fermi, who used a purely classical definition of the angular momentum, 


L=rxp (4.1) 


and applied a thermodynamical treatment imposing, besides energy conservation, the 
condition that the total angular momentum should have a prescribed value. For this 
approximation to be valid, it would be necessary to go up to a very high energy and to have 
a large number of secondary particles in equilibrium. Even then it is not clear if, or under 
what conditions the quantum effects of angular momenta would turn out to be really negligible. 

In our treatment the quantum-mechanical counterpart of relation (4.1) is found in the 
cut-off factor (prj(pr))? in the phase-space integrals. Consequently, our method can in 
principle be applied to more detailed statistical treatment at lower energy. It is expected that 
the application of our formalism to higher energies will be able to clarify the question of 
the validity region of the thermodynamical approximation. 

Further, we do not regard the so-called Lorentz-contraction of the effective interaction 
volume in a literal way. This deformed shape of the volume appears, on the other hand, 
to play an essential role in Fermi’s treatment. 

4. Another difference between Fermi’s approach and ours is that in the former the 
inelasticity of the collision is always put equal to one and the angular momentum of the 
system is specified by the classical impact parameter of the collision, while in the latter 
the total energy W and the total angular momentum L, M, of the assembly of secondary 
particles are regarded as certain independent parameters, without entering into the 
mechanism how they have been realized. 

From the experimental data of super-high energy collisions (cosmic ray jets) we have 
the feeling that at least two independent parameters would have to be introduced in order 
to classify the events іп a sensible way!9. As these two parameters we take here the most funda- 
mental quantities: energy and angular momentum, But the purpose of the proposed statistical 


19) In a definite energy region, for example, we find quite a variety of events (isotropic distribution of 
large number of secondary particles, very anisotropic distribution of small number of particles, and so on) 


among quasi-nucleon-nucleon collisions (those collisions which appear essentially similar to nucleon-nucleon 
collisions). j 
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analysis is, for the time being, not to answer the question: “How does this transfer of the 
energy and the angular momentum from the incident to the secondary particles take рїасе?”, 
but only: **What kind of number and angular distributions are expected, if we assume this 
transfer of energy and angular momentum 2”? 

_ 5. А sharp criticism has been raised against Fermi's application of the statistical theory 
to long-distance collisions (Feinberg and Chernavski 1951; Fukuda 1953; Rozental and 
Chernavskii 1954). The point is that for a distant collision the statistical equilibrium could 
not be established without admitting a propagation of disturbances with super-light velocity 
because the two particles fly past one another in a very short time. Thus the problem here 
is concerned with the mechanism of the peripheral collision. 

We have already emphasized, however, that we are not really discussing the mechanism 
of the collision. In other words, we try to separate from the general expression for the transi- 
tion probability 

transition probability ~ [matrix element|? x (phase space volume per unit energy), 


the second factor, which represents so to speak the “kinematical? aspects of the 
collision. In order to investigate the latter effects we assume the simplest mechanism: the 
square of the matrix element is proportional to the probability that the final free particles 
are found simultaneously in the effective volume of interaction. 
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M. Danysz, Professor M. Miesowicz, and Professor J. Werle for their hospitality extended 
to me during my stay in Warsaw. Further, I am grateful to many of my colleagues in Warsaw, 
especially to Dr P. Zieliñski and Mr. R. Raczka, for their comments, and also to my wife 
Sachie Koba, who helped me in computation. 


APPENDIX 1. 


Approximate evaluation of J). 
We have to evaluate expression (3.3): 


оо 


225 1 н z 
ج‎ d R—— w % . Al 
1-25 1 2) екр(— iau) (АЛ) 
IR 
(i) For 1 — 0, 
h= ах dp exp (— ia Vp? + n°). (A.2) 
0 


Putting p = ш sinh z, w = Vp? +u? = p cosh 2, we get 


h= 2g Ен E cosh 2 exp (— iam cosh t) = — 69 Ru Hf? (au) (А.З) 
л 


0 
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If au < 1, which corresponds to the extremely relativistic case, 


ДРОМ. 
H9 (ар) = 2. — 


аш” 
so that 
ее, (А.4) 
” 1% 


The same result may be obtained directly from (A.2) by replacing w by p. 
If on the other hand «u2»1 (non-relativistic approximation), we have 


1 H (au) ~ = exp (— i&u) exp (- 1 i) , (А.5) 


ig Vs exp | —i z) гек“ exp (— iam). (А.б) 


The same result may be obtained directly from (A.2) by replacing w by u+ p?/2u. 
(ii) For 12» 1 the integration domain of the variable p is always greater than 1/R = w 


whence 


so that we can apply the extremely relativistic approximation. 
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Here we have used the relation au«£1 and made the substitution p —l[Rzq. 
To carry out the integration in a closed form, we make a further approximation by 
replacing the function (I/R)(q-4- 1|R)-* by another function of similar behaviour”; 
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20) Both of these functions have the value опе ай 


q = 0 and decreasing monotonously with g, they tend 
to zero with q — co 
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APPENDIX 2. 


Evaluation of 5,11, Jy, ..., Ly). 

We insert the expressions for 1, which have been obtained in Appendix 1, into 
(3.2) and carry out the final integration over the variable g. 

(I) When all the particles are treated as extremely relativistic, the problem is reduced 
to partial fraction decomposition and evaluation of the residues, and can be worked out 
in an elementary way. For example: 
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(П) When an even number of non-relativistic S-wave particles are participating in the 


equilibrium, the situation is quite similar to the above case. For example: 
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(III) When ап odd number of non-relativistic S-wave particles and no other particles 
of highe- angular momentum are present, the problem may be again simplified as follows: 
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(IV) When we have an odd number of non-relativistic S-wave particles as well as some 
other particles of higher angular momenta, we may proceed as shown by the following 


example 
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Here the integrand s,(...) means the phase space volume for three particles with available 
energy W — u — p*/2u instead of W, and is regarded as a function of p. Thus 
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where F is defined by 


x 
F(x)-— en fier" dy. 
2 0 
When we have to do with non-relativistic particles, the argument of F ік relatively small, 
so that the numerical value of F can be easily obtained. 


APPENDIX 3. 


Table of relative magnitude of phase space volumes for the case RN = 10 u, 9 = 7/8 
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Appendix 3. (Continued) 
L= 0 L-1L22 Г--8 L=4 L=5 L6 L=7 L=8 1-9 


(2,5) шэг e 7 7 7 2 7 = 
(2,6) хэн = == == 2 2 2 2 2 — 
(2,7) s КЕЗЕ m — = 0 0 0 0 0 
(3,3) 8 8 8 8 8 8 E 
(3,4) 6 6 6 6 6 6 6 — MM 
(3,5) әнді тай 2 2 9 2 д > омж 
(3,6) — — -- 0 0 0 0 0 0 0 
(4,4) 1 1 1 1 1 1 1 1 n= 
(4,5) = 0 0 0 0 0 0 0 0 0 
(0,0,0) 52 - — 
(0,0,1) — 709 — — 
(0.0.2) L3 чүр — = 
(0.0.3) "ырь ا کے‎ “бар - = 
(0.0.4) 20 

(0,0,5) EMEN UI DO 
(0,0,6) = ар l — — 
(0.0,7) — e F УЕ o 
(0,1,1) 59 59 59 -- “Қ 
(0,1,2) не 1727; 04: 4401-2125 c о 
(0,1,3) —— 2 хээл rag 159: 99 = 
(0,1,4) 197 £9 Жо 

(0,1,5) = 5 5 5. = АН 
(0,1,6) یر و‎ ۰ D ИЛ 1 1 - 54 
(0,1,7) 0 0 0. Ж 
(0,2,2) 15. IS 46 515 5 — RN 
(0,2,3) = IB. 48 Шз Б а 2 
(0,2,4) Ap 4 4 4 4 4 ---- _ 
(0.2,5) e ا‎ 1 1 1 1. У 
(0,2,6) Қы ey C MEE 0 0 0 0 0 = 
(0,3,3) 2 2 2 2 2 2 2 ی‎ 
(0,3,4) т 40 мй 10 Q0 02 аа 
(0,3,5) vr с. 0 0 0 0 0 0 0-2 
(0,4,4) 0 0 0 0 0 0 0 0: = 
(1,1,1) 10: 30 Юзи ани s ЫИС 
(1,1,2) 13 26. мо "006 16 E гүтгэн 
(1,1,3) — 5. 410, 91652810 Т 
(1,1,4) == — 2 4 6 4 2 == e = 28 


(1,1,5) 
(1,1,6) 
(1,7) 
(1,2,2) 
(1,2,3) 
(1,2,4) 
(1,2,5) 
(1,2,6) 
(1,3,3) 
(1,3,4) 
(1,3,5) 
(1,4,4) 
(2,2,2) 
(2,2,3) 
(2,2,4) 
(2,2,5) 
(2,3,3) . 
(2,3,4) 
(3,3,3) 


(0,0,0,0) 
(0.0,0,1) 
(0,0,0,2) 
(0,0,0,3) 
(0,0,0,4) 
(0,0,0,5) 
(0,0,0,6) 
(0,0,1,1) 
(0,0,1,2) 
(0,0,1,3) 
(0,0,1,4) 
(0,0,1,5) 
(0,0,1,6) 
(0,0,2,2) 
(0,0,2,3) 
(0,0,2,4) 
(0,0,2,5) 
(0.0,3,3) 


Appendix 3. (Continued) 


221 


OLN І--2 ІГ--377,-4/ L5 1-6 L=7 L= 8 L=9 


== — -- 0 1 || Jl 0 — — 
lI — == — 0 0 0 0 0 — 
БЭР 22 = == == 0 0 0 0 0 
5 14 14 14 10 5 — — — — 
3 6 9 9 9 6 3 
= 1 1 2 2 2 1 1 == — 
== == 0 0 0 0 0 0 0 — 
18 == Зэн 0 0 0 0 0 0 0 
0 1 1 1 1 1 1 0 -- == 
0 0 0 0 0 0 0 0 0 — 
-5 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 1 2 2 1 1 0 
0 1 1 1 1 1 i 0 == — 
0 0 0 0 0 0 0 0 0 LÍ 
-- 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 — 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
ТЛ! 
e 23 = ET 
Ж. кз 12 
3 
hk 2 -- E 
0 
0 
16 16 16 Е 
== 14 14 14 E 
‚Ж, x. 4 4 4 — -- = 552 БЕ 
1 1 1 = 
27 #2 Et = 0 0 0 
E i 0 0 0 -- = 
2) % 25-22 2 
27 1! 1 1 1 1 
к = 0 0 0 0 0 
0 0 0 0 0 == E 
0 0 0 0 0 0 0 


232 


(0,0,3,4) 
(0,1,1,1) 
(0,1,1,2) 
(0,1,1,3) 
(0,1,1,4) 
(0,1,1,5) 
(0,1,1,6) 
(0,1,2,2) 
(0,1,2,3) 
(0,1,2,4) 
(0,1,2,5) 
(0,1,3,3) 
(0,1,3,4) 
(0,2,2,2) 
(0,2,2,3) 
(0,2,2,4) 
(0,2,3,3) 
{1;1,1,1) 
(1,1,1,2) 
(1,1,1,3) 
(1,1,1,4) 
(1,1,1,5) 
(1,1,2,2) 
(1,1,2,3) 
(1,1,2,4) 
(1,1,3,3) 
(1,1,3,4) 
(1,2,2,2) 
(1,2,2,3) 


Appendix 3. (Continued) 


L=0 L=1 L22 L3 L=4 L=5 L-—6 L= 


— 0 0 0 
1 12 8 4 
4 9 13 9 

-- 2 3 

-- - 0 1 

0 
1 3 3 3 
0 1 1 1 

-- 0 0 0 

— -- 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
1 2 2 1 
1 2 2 2 
Ü 0 0 0 

= 0 0 0 

— — 0 0 
0 1 1 1 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
1 — = = 

= 9 ҮР ТҰ 

> == ih p 


L8 Lem 

0 0 0 0 = ш 
4 25 
2 ii — == === =. 
1 1 0 

0 0 0 0 = = 
0 0 0 0 0 — 
2 1 A 
l 1 0 —  — 
0 0 0 0 == = 
0 0 0 0 0 = 
0 0 0 0 — I 
0 0 0 0 0 -- 
0 0 0 => = цэн 
0 0 0 0 = Жл 
0 0 0 0 0 — 
0 0 0 0 0 е 
0 a 
1 0 2-- 

0 0 0 = = 28 
0 0 0 0 хээ = 
0 0 0 0 0 — 
0 0 0 

0 0 0 0 a AR 
0 0 0 0 0 Et: 
0 0 0 0 0 = 
0 0 0 0 0 0 
0 0 0 0 — 5 
0 0 0 0 0 E 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 кы 
0 0 0 0 0 0 


is 
© IO! 
=з SIS! 
© S! ol 


5 
m 
D 
m 


m 
EJ 


B 
5 
m 


5 
` 


Si o> CO! 


2) 


= = еш 


heu ele» 


2 


Da r ف‎ ык St UA 


СООТ СЭТ Ee» Ss 
Һ9 59 fem - ык = = 


5 
m 
э 


5 


ет: 
© OI! O 


e 
ee 
со 
— 


EERS 


5 


(0,0,1,1,4) 
(0.0,1,1,5) 
(0,0,1,2,2) 
(0.0,1,2,3) 
(0,0,1,2,4) 
(0,0,2,2,2) 
(0,0,2,2,3) 
(0,1,1,1,1) 
01112) 
(0,1,1,1,3) 
(0,1,1,1,4) 
(0,1,1,2,2) 
(0,1,1,2,3) 
(0,1,1,2,4) 
(0,1,1,3,3) 
(0,1,2,2,2) 
(0,1,2,2,3) 
(0,2,2,2,2) 
11:1,1,1,1) 


Appendix 3. (Continued) 


233 


=0 L=1 L=2 L=3 L4 L=5 L=6 L=7 L8 L=9 


OS coc c с Sb. 


wc Ба СЭ | 


Soo со же с с с есе гс = © || 


Se © (е) (зэ) (е) бэ? (22) Kay (аз) (е) (е) Se (22) aa (=; 


о С н 


| 


0 = ER 
2 0 = m Pa 
0 Я 

0 0 

0 0 0 

0 S 

0 0 

0 0 0 

0 0 0 

0 

0 0 = 3 xus) (б—— 
0 0 0 

0 0 0 Оре 
0 (055 15/558, 5с = ra eS 
0 0 OSA 
0 0 0 Q тн 
0 0 oe - 
0 0 0 Са == 
0 2223 == 

0 0 m 
0 0 0 cA 1-2 Жм: 
ООУ Krone ДИ oie Me: 
0 0 (уеде = лу алдлан 
0 0 0 О. Hr 
0 0 0 0 Өл Ж ор 
0 0 0 0 0 = 
0 0 0 ы ыы ج‎ 
0 0 0 0 О 
0 0 0 0 0 = 
0 Охна 


234 


Appendix 3. (Continued) 
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Configurations with N > б are not listed here, because they do not give appreciable 
contributions. The remarks to Table 1 apply also to this appendix. 


NOTE ADDED IN PROOF. 


1) Aformulation similar to this work has been presented by Cook and Lepore in connection with the 
analysis of nucleon-antinucleon annihilation. See Cook L. F. Jr. and Lepore J. V., Phys. Rev., 120, 1028 
(1960). Evaluation of the number of states with angular momentum conservation (but without energy 
conservation) was once discussed by Magalinskii. See Magalinskii V. B., Zh. eksper. teor. Fiz., 36, 93 (1959). 
I am indebted to Dr. V. S. Barashenkov for calling my attention to the latter work. 

2) General properties and methods of evaluation of the weight factor Z(L, Û, ..., ly) which appear in (2.17) 
and (2.18) have been discussed elsewhere. See Koba Z., “On the angular momentum weight factor in the statistical 
theory of multiple particle production. 1. and IL”, Bulletin of the Polish Academy of Sciences, in press. 

3) In (2.10) and the following expressions it is tacitly understood that the summation over 1, 8 are taken 
under the restriction | < l, <... < ly. See the combinations appearing in Table 1 and Appendix 3. Hence we 
have the factor (лу! n, ! ... n, !)-! in order to take into account the indistinguishability of pions. For later pur- 
poses it is more convenient, however, to remove this restriction and take all possible values of 1,78. Then we have 
simply the factor J/N!, just as in the case of plane wave representation. In derving (2.21) this relation has been 
used. 

4) The short summary of the present work, which appeared in Nuovo Cimento, 18, 608 (1960), contains 
unnecessary double summations in (13) and (14), the counterparts of (2.20) — (2.22) of the present text, for 


angular distribution of secondary particles. Those double summations should be reduced to single ones by putting 
т=т/. 
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ON THE BREADTHS OF ANNIHILATION LINES IN ONE- 
AND TWO- VALUED METALS 


By BRONISŁAW SREDNIAWA 
Institute of Theoretical Physics, Jagellonian University, Cracow 
( Received September 29, 1960) 


Stewart's measurements (Canad. J. Phys., 35, 168 (1957)) of the breadths of annihilation 
lines in the process of annihilation of positrons in metals are analysed using the method of thet 
author (Helv. phys. Acta, 33, 131(1960)) in the case of 16 one- and two-valeud metals. The resuls, 
are: for light metals (Li, Be, Na, Mg) positrons are annihilated with conduction electrons only 
for heavier metals the contribution of annihilations with the electrons of the atomic shells in the 
ions increases gradually. Noble metals have a particular position, for Cu, Ag, and Au the annihila- 
tion takes place almost exclusively in atomic shells. 


1. Introduction 


In the paper entitled “On the breadths of annihilation lines in copper and in gold” 
(Sredniawa 1960, henceforth referred to as 11), a certain method of theoretical investigation of 
the breadths of annihilation lines of y-rays arising in the process of annihilation of electrons and 
positrons in metals has been proposed. The positrons produced in the f-decay of Сиё (with the 
maximal initial energy of 0.66 MeV) are slowed down in a sample of metal and then annihilate 
with electrons. The process of slowing down may end 1) outside the ions of metal lattice or 
2) inside the ions. In the first case the velocity of the positrons falls down to thermal velocities, 
which we neglect, assuming that the positron is at rest. This positron annihilates then with 
a conduction electron. In the second case the positron is stopped inside the ion and annihi- 
lates with one of the electrons of atomic shells of the ion. 

In I theoretical line breadths of the annihilations in the regions of different electronic 
shells of the atom of the annihilation with conduction electrons are calculated. Contributions 
from the process of different regions where annihilation may take place to the real annihilation 
process аге then chosen in the manner to get by their superposition such a shape of the 
theoretical annihilation line which would be the best approximation of the experimental 
results. Since wearetempting to examine the consequences of this assumption on the possibility 
of annihilation of the positron inside the ion we do not consider the second alternative, 
namely the excluded volume effect (de Benedetti et al. 1950, Ferrel 1956) based upon the 


assumption that the positrons cannot enter inside the ions. We neglect the possibility of 


positronium production and the effects due to plasma, which according Ferrel do not play 


a great role in this process. 


1 We accept notations of this paper. 
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This method was applied in I to copper (comparison with experimental results of du 
Mond et al. (1949)) and to gold (comparison with experimental results of de Benedetti et al.). 
'The present paper presents a more detailed analysis of the results of the measurements of 
Stewart (1957). These measurements performed by means of the method of angular 
correlations of the photons produced in the process of annihilations comprise the greatest 
number of elements. Our analysis concerns one- and two- valued metals for which (in 
a certain approximation) the conduction electrons may be regarded as free and represented 
as plane waves. We have analysed 16 metals of groups Ia, Ib, Па, ПЫ, of elements. 
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Fig. 1 
—— experimental curve 
EE эе theoretical curve, L shell 
---- “-- theor. curve for cond. electrons 
ОЛОРОЛЛОР resultant theoretical curve (in tables 2, 3, 4)? 


3 For table 1 see discussion 


237 


2. Analysts 


The starting point for the analysis are the experimental coincidence counting - rate 
diagrams of Stewart’s paper (p. 171). From this diagram, using the method of de Benedetti 
et al. the average momentum p of the mass center of annihilating particles and then the 
reduced average momentum of the mass center D = p[mc was determined. Stewart's 
measurements were not corrected for the finite length of detector slits, but de Benedetti's 
method was applied up to abscissae of the curves of Stewart not greater than 5 miliradians 


where this correction does not play a significant role. Having Р, 


formula (1, 85, p. 139) 


с. We calculate, using the 


P... = 0.082 о (1) 


the dispersion o of the gaussian curve corresponding to this average momentum. Results 
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of these calculations for 16 analysed metals are given in columns 3 and 4 of table 1. “Ехреті- 
mental” curves are drawn in diagrams 1—4. These curves, like all other on these diagrams 
are normalized to the unit height. The maximum energy ву of the conduction electrons 
(of Fermi gas in absolute zero temperature) are taken from Seitz (table XXVIII, p. 146 and 
formula (20), р. 145). They are given in table 1, col. 5. Having & for considered metals 
we calculate corresponding reduced momenta P' (see I, $ 4) and draw the curves (normalized 
to height 1) of annihilation of conduction electrons with positrons at rest, according to 
formula for I (1, p. 146), where instead of K we put as variable 


3 

(norm) | |44) 

13 -1- (aa Va (2) 
“0 

Іп order to get the contribution of annihilations іп ions we take the maxima of the charge 


density from llartree diagrams of tables (references are given in Gombas 1950 and 1), from 
which according to (I (21) and I (26)) we calculated the annihilation curves for different 
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Table 1 
- ч | 
Metal 2 Уд 1020 £g (eV) Metal 2 LOPATA 1020 eg (eV) 
Li 3 2,6 3.2 4.72 Rb 37 8.2 3.9 1.78 
Ве 4 4.0 4.9 14.3 Sr 38 4.6 5.6 3.85 
Na 11 2.2 DT 3.12 Ag 47 7.4 9.0 5.51 
Mg 18 8.6 4.4 ТЛА Са 48 4.6 5.6 7.53 
K 19 Qu 3.3 ОЛА Сз 55 4.1 4.1 1.54 
Ca 20 3.4 4.1 4.26 Ba 56 4.9 4.9 3.65 
Cu 29 8.0 9.7 7.04 Au 79 7.6 7.6 5.54 
Zn 30 4.2 5:2 9.60 Не 80 8.0 8.0 6.80 
l 
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shels. For Cd, Ba, Hg the same maxima of charge density were assumed as for their neigh- 
bours in the table of elements Ag, Cs, Au respectively. 

Then by trials a curve (superposition of the above mentioned theoretical curves) was 
chosen which was the best fit of the experimental line. In table 2 the contributions of different 


electron groups in the process of annihilation are given. 


Table 2 


Results 
Contribution of (in %) 


Group Metal Shell L | Sh. М Sh. N Sh. O Cond. el. 
Ia Li 100 
Na 100 
K 30 70 
Rb 50 50 
Cs 20 80 0 
Ib Cu 5 95 0 
Ag 20 80 0 
Au 20 80 0 
Ila Be 100 
Mg 100 
Ca 40 60 
Sr 10 90 0 
Ba 40 60 0 
IIb Zn 10 90 
Cd 70 30 
Hg 30 70 0 


3. Discussion 


Because of crude approximations assumed in I, the contributions given in table 2 have, 
of course only orientational meaning. We see that for light metals Li, Be, Na, Mg, (fig. 1) 
positrons annihilate almost exclusively with conduction electrons. Therefore the “resultant 
theoretical curves" are for them identical with “experimental curves”. Let us now cut 
the tails of experimental curves for these metals in order to imitate the theoretical curves 
of Fermi distribution. Taking the abscissa of the point, where these curves cut the A/-axis 
and taking account of (2) we get a crude experimental value of maximal Fermi energy 
for these metals (see dotted lines in fig. 1): 


Table 3 


Metal | Li | Ве | Ха | Mg 
(еі) еке (e V) | 3.5 | 11 | 1.6 | 74 
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These values for Li and Mg approach to values calculated from energy bands (see Seitz 
1940, p. 440), for Be and Na they are much smaller. For heavier metals the contribution 
of electrons of the external atomic shells increases, more quickly for metals of the group Ha 
than for their neighbours of the group Ia. 

Particular position in this table have metals (group I5) Cu, Ag, Au, where (according 
to our assumptions, given in the introduction) the annihilation takes place in exterior shells 
of the ions. In the neighbouring metals of the IIb group the contribution of conduction 
electrons is much more greater, except of Hg, which behaves similar to Au. 

Let us orientatively take for the radius of an atom the distance in which the charge 
density falls down to the l/e-value of the most external maximum. Let us calculate the 
volume occupied by ions in 1 em? of the metal sample and assume that the probability of 
annihilation of the positron by a conduction electron (or electron in an ion) is proportional 
to the volume of empty space in the metal sample (or volume occupied by ions) and the 
density of conduction electrons (or density of electrons in external shell). Then we obtain 
the ratio JV of probabilities for annihilation in ions and outside ions given in table 4. 


Table 4 
Metal | Cu | Zn | Ag | Cd | Au | Hg 
W | 2.5 | ыа 11 | 4 | 5 mee: 


In these calculations we did not take into account the repulsion of positrons by electro- 
static fields of nuclei (almost equal for neighbouring metals), nevertheless probabilities 
of annihilation can be indeed much greater for noble metals than for their neighbours. 
The results for noble metals and heavier metals of group IIb may also be interpreted 
in such a way that the simple Sommerfeld model of the metal consisting of ions and conduc- 
tion electrons forming a Fermi gas is not valid and the conduction electrons have energies 
not very different from the energies of bound electrons in outer shells of ions. In this case 
we cannot distinguish distinctly between these classes of electrons and their contributions 


to the annihilation process. 
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The real and imaginary part of the medium and high energy potential for nucleon-nucleus 
scattering has been calculated for an infinite nuclear medium with help of the recent phenom- 
enological nucleon-nucleon phase shifts. For the S wave contribution the exact solution of the 
equation for the nucleon-nucleon scattering for the separable nucleon-nucleon interaction has 
been used. In calculation of the I > 0 contribution to the imaginary part of the optical potential 
the Goldberger method was applied. The results are compared with experiment and discussed. 


$ 1. Introduction 


The present paper deals with the derivation of the parameters of the optical potential 
for the nucleon-nucleus scattering based on our knowledge of the nuclear forces. 

Тһе general theory of the optical potential has been given by Watson and his collaborators 
(Watson 1953, Francis and Watson 1953, Takeda and Watson 1955). In the last few 
years several calculations of the optical model parameters have been presented; all of them 
have been based on Watson’s theory. According to this theory the optical potential V for 
a nucleus with high mass number 4 is given by the following expression: 


A 
Бы А (1) 


where yy is the ground state wave function of the target nucleus. The transistion operator to; 
describing the scattering between the particle ,,0” and „й? in the presence of the other 


target nucleons satisfies the integral equation: 


toi = Voi + Voi = Loi; (2) 


(245) 
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where the operator Q excludes all the states occupied by the target nucleons, and 


е= b= Hi li Ad 2) (3) 


where vg; is the nucleon-nucleon interaction, E is the ground state energy of the target 
nucleus plus the energy of the projectile, H, is the hamiltonian of the target nucleus and 7, 
the kinetic energy of the projectile. 

The expression (1) for the optical potential V already contains some approximations: 

a) It is understood that one has to antisymmetrize the 1; matrix with respect to the 
particles o and 1. This takes into account the “primary exchange effect" in the terminology 
of Takeda and Watson (1955). The so-called “target exchange effects" are neglected in this 
way. Тһе estimation of these effects by Takeda and Watson (1955) [compare also Sawicki 
(1960)] indicate that for energies high enough (more than 50— 100 MeV) these “target 
exchange effects" do not play any essential role. 

b) In expression (1) the incoherent scattering is neglected. This incoherent scattering 
representing subsequent inelastic scattering of the projectile by a few different nucleons 
is difficult to estimate. The estimates by Watson (1953) and by Dabrowski and Sawicki 
(1959) indicate that the incoherent scattering correction smay be neglected for about 
the same energies for which the limited antisymmetrization described in a) is justified. 

Both simplifications a) and b) restrict the validity of the expression (1) for medium 
and high energies (z 50 — 100 MeV). 

In applying eqs. (1) — (3) of Watson's theory in actual calculations of V one is, however, 
forced to introduce further simplifications. 

As a rule one neglects the imaginary part of V in eq. (2) for the ¢ matrix, i.e. one puts 
V лг Re V in eq. (3). This approximation is hard to justify in a frame of a theory valid for 
high energies where the energy dependent imaginary part of V, Im I, is of the same order 
of magnitude as Re / [see, e.g., Fig. 9 in the review article by Glassgold (1959)]. 

The energy dependence of Re V in eq. (2) for the 4 matrix is usually approximated 
by the quadratic dependence of the momentum. This is equivalent to the replacement of 
the nucleon mass M by an effective nucleon mass M*. Even if one goes beyond the quadratic 
approximation one still can use the effective mass, which depends then on the momentum. 
(Dabrowski, Sawicki 1961). 

With all these approximations it is still difficult to solve the equation (2)for the £ matrix 
in the case of realistic nuclear forces. Verlet and Gavoret (1958) make use of a separable 
potential which fits the scattering data of low energies (Yamaguchi 1954). With this potential 
it is easy to solve the ¿matrix equation. (2). However, the Yamaguchi potential acts 
in the 5 state only and its application in the medium and high energy optical model caleulation 
may be questioned. Namely, at higher energies the higher partial waves play an essential 
role in the nucleon-nucleon collisions. On the other hand, the Pauli principle, i.e. the Q 
operator in eq (2), is less important at higher energies. Namely, at these energies the inter- 
mediate states excluded from the ¢ matrix equation play a relatively smaller role. This fact 
may be inferred easilly from the numerical results obtained by Verlet and Gavoret (1958— 
their Table 1) and by Dabrowski and Sawicki (1960 — their Table ПІ) where the optical 
potential without the Pauli principle (Q — 1) is also considered. 
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Кот the reasons mentioned above, we feel that опе improves the calculation of the 
optical model parameters by including the higher partial waves into the £ matrix neglecting 
completely the influence of the Pauli principle on these higher partial waves. This program 
may be easily realized as the t-matrix satisfying now eq (2) with the operator О = 1 may 
be expressed with the help of the nucleon-nucleon phase shifts. A calculation of this 
type of the imaginary part of the optical potential has been presented recently by Gomes 
(1959). 

In the present paper we replace V in the denominator of the £ matrix equation (2) by 
an energy independent real potential. Thus we restrict ourselves to the case M* — M 
because only in this case by dropping the Q operator in eq. (2) does one obtain the equation 
for the scattering of two free nucleons. 

Since the Coulomb effect is not taken into account in this paper the results obtained 
are aplicable for an incoming neutron rather than a proton except at the highest energies 
considered. 

In the next section (8 2) the general expressions for V in terms of the nucleon-nucleon 
phase shifts are given. In the calculation of ImV we follow the Goldberger method known 
also as “frivolous model" (Lane, Wandel 1955). Thus we include some of the effects of 
the exclusion principle. One can show that, this, ‘frivolous model", or — in the terminology 
of Gomes (1959) the “final states” approach, may be obtained as an approximation of Watson's 
theory (Dabrowski, Sawicki 1961). 

In 83 the contribution to V from all the partial waves with the orbital quantum numbers 
1 > 0 are calculated and added to the results of the S wave t matrix of Verlet and Gavoret 
(1958). We use the phase shifts corresponding to the Gammel, Christian and Thaler (1957), 
Gammel and Thaler (1957), and Signell and Marshak (1958) potential. The results are 


then compared with the experiment and discusssed. 


8 2. General expressions for the optical potential in terms of the nucleon-nucleon phase 


shifts 


In the following we restrict ourselves to the Fermi model, i.e. to the nuclear matter. 
Namely, we assume that the ground state wave function yy of the target nucleus is the anti- 
symmetrized product of plane waves. It is easy then to calculate in momentum representa- 


3 
CV] К) = Ё 212222 Xs (To) X 


xe (с) y, (1o) t Te (4) 
where 7), y, are the spin and isospin functions; s and t are the corresponding “magnetic 
numbers". With the assumed form of y, one gets easily (Dabrowski and Sawicki 1959): 


«ЕК? = V(ko) (ко — ko) (5) 


V (ko) = 4 Ї dkt 2-5 (6) 


226432 


Bon, i.e. 
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where i(k) = «ККУ is the on — the — energy — shell forward scattering component of the 
average t matrix. The average t matrix appears because of four nucleons (neutron and proton, 
each of them with spin up and down) occupying each state of a given momentum k. In the 


following it is assumed that the incoming nucleon is a neutron. For the ¢ matrix one gets: 


«ККУ = MK 28192) (7) 


ЖБЕК К 


where 
CK'|t,,| E> = [CR |1- 
t_ ko CRE 


3t— 


ЭЕ 5 
ЭЖ (8) 


where the indices 1, 3 апа +, — denote the scattering in the singlet, triplet spin, and in the 


КУ], 


k >t: ОК' 


1+ 


CK 


k»,-- 2 Ck! | 


even, odd orbital states, respectively. Furtheron we have for each of the 535: 
k»,—$[CK'|t k>] 
К>, =4[<k'|t|k>+<¢—k'|t|k>] (9) 


According to Watson's theory each of the 53 ¿£ matrices is supposed to satisfy the equation 
(2) with the corresponding 135, 


CR t 


ку - 4-Е 


«к 


An exceptional case for which one can solve eq (2) for the t matrix easily is the case of 
the separable nucleon nucleon potential: 


1:43 d жый 
(as I "арте" = «роз 


p> 


=~ ep!) в(р) (10) 


Depending on the form of the g function we call the potential (10) the “Yamaguchi poten- 
tial" (Y)(Yamaguchi 1954) when 


8(р) = (p?+ B?)-1 (11) 


and the “effective range potential? (E R) (De Dominicis and Martin 1957, Martin and 
Gourdin 1957) when 


g(p) = (p? + В2)-% (12) 


The parameters 2, В have in the singlet and triplet states the respective values 14, 18:31:38 
adjusted to the low energy nucleon-nucleon scattering data. 

Although the separable potential (10) allows us to solve the £ matrix equation easily, 
it has disadvantages; first of all this potential acts in the S state only. According to our 
remarks in the introduction we want to improve the results for the optical potential obtained 
with the separable potential by adding the contribution from higher (than 5) partial 
wawes in the phase shift approximation. To get the phase shift approximation for the 
t matrix we put the operator Q = 1 in eq. (2) and replace V with an energy -independ- 
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ent real potential. In this way our £ matrices go over into the ty free nucleon transition ma- 
trices, i. e. putting Q = 1 in eq (2), is equivalent to the approximation: 


adi (13) 
It is well known how to express the on — energy — shell elements of the tọ matrix with 


help of the nucleon-nucleon phase shifts. The relevant expressions are given, e. g., by 
Stapp et al. (1957)1, As the result we get: 


to(p) = Re to(p) + Im to(p) (14) 
where 
Re 50) NES n. (15) 
Im to (p) 32z2 Мр \2/(р) 
where 


g(p) = [ >) (oddj) + 3 >) (even j)] (2j t 1) sin 20; 4- [У (even j) + 3Х(о44/)| х 
x [(Zj + 1) sin 25; + (2j + 3) (sin 25; |, a + sin 29, , „)] + 3 sin 27 


where 6; denotes the singlet state phase shift and 7j, 7j; a> 7j; в denote the uncoupled (/ = j) 
and respectively coupled (J = j — 1, / = j + 1) triplet phase shifts of Blatt and Bieden- 
harn (1952а, b); 7, is the 3P, phase shift. The equation for /(р) is obtained from eq. (16) 
by replacing sin 20, sin 2) by sinê, sin?r, respectively. The above equations for ty may be 
found in the papers by Riesenfeld and Watson (1956) and Ohnuma (1958). The phase 
shifts are usually given as function of the laboratory energy E, related to p by E, = 2h?p?/M. 

To calculate the optical potential V(kj), eq. (6), let us introduce in analogy to Verlet 
and Gavoret (1958) the coordinates: 


In the new coordinates we have 
Got kee (сэрэг 
dk, = чя [ duu (19) 
Ер бар) max Гр), үйе 
and for V(ko) we get 
(y+ ЕЕ)? V disks p 


Re V/(kj) 2h? &(p) | 
a vi 77 za, | Marl eem Ей, 


(ko—kp)/2 max (а.-», к-р] 


Ё l1 
2n* M 2 
forward scattering matrix: M, = 3 (2My + Moo M Janen 18 Table III of Stapp et al. (1957) one 
finds the expressions for M,,, Moos Mss , for the п-п case by dropping all the Coulomb factors. In the n-p case 
the expression for М1, Moo, Mss may be obtained from the corresponding expressions for the n-n case by divid- 
ing them by 2 and by dropping all restriction imposed on the / summations. 


1 To obtain our eqs. (15), (16) one has to notice that our E. [Mint M,ple=o> where the 
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If we introduce the dimensionless variables v = kọ/kp, x = p/kp we get finaly: 


(T 1)/2 


Ке V(k)\ 1 2 aka 1 É " (-- 211 НУ | | 
Im Fik) xz M v ede 2 2f (x) 9% 


(v—1)/2 


Тһе eq. (20) is valid in the case k; > V2 kp, for which (k, — p) > И, — p? — p? and for which. 
we can replace the lower limit of the u-integration by (k; — p). The connection between: 
Ej and x is now: 

hk m 


Ет, -- 2 М 


Now we want to improve the expression for Im V (kọ) by taking into account some of 
the effects of the exclusion principle. For this reason we follow the Goldberger (1948) 
method, known also as the “frivolous model" (Lane and Wandel 1955; Clementel and. 
Villy 1955). The relation of this method (and its generalization for an arbitrary effective. 
mass) to Watson's theory of the optical model is presented in the paper by Dabrowski and 
Sawicki (1961). Here we follow a more intuitive approach. 
2 
| 3 
ау 


We start from the optical theorem: 
where the subscript “ау” denotes the proper average over spin and isotopic spin states, 


~ aM 
Im (р) =— СА) 40, Ke пыр) 


According to the Goldberger method we exclude from the Q, — integration those 
parts that would lead to final state (q/q)p occupied by other nucleons in the Fermi sea. 
For this reason we introduce the О operator defined explicitly by: 


1 q p|- 1 for 
0|2 + - Ky) + 2 


О otherwise 


1 q 
> (ko + Ку) + PET (22) 


and replace Im te p) by Im EEE (p): 


- zM 
Im £& (+k) (p) = — on: Р Гао, x 


1 
x Q Ë (kot Kj) + 1| Ke 72) 


If we neglect the angular dependence of ЇЕМ) 


1 e 


2 
| ‚ we may rewrite eq. (23), 
av 


using once more eq. (21), in the form: 


Im nx) UD) = Im tolp) лс _ fao, х 


T E (otk) 3 1 ян” 
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Since ||? is proportional to the differential cross section, we have neglected angular 
dependence of the differential cross section while going from eq. (23) to eq. (24). Stricly 
speaking we have neglected the terms proportional to P, (cos Ө), P,(cos O), ... (O is the 
scattering angle) in the expansion 46/40 = XA,P, because the Q integration leaves 
the even / terms only. The P, term was evaluated by Verlet and Gavoret (1958), who found 
that it results an increase of the Im V of an order of magnitude of 10 per cent. At first 
sight the smallness of this correction is surprising in view of the anisotropy of the neutron- 
proton scattering. However, we must not forget that the Q operator in eq. (23) excludes 
scattering with a small momentum transfer and in the £2, integration region, which is 
not affected by the Q operator, the anisotropy of the neutron-proton differential cross 
section is less pronounced. 

In the variables u, p, eqs. (17), (18). we obtain easily: 


(hot kp)? | ik» 


Ол" 
Im Voa(ko) = — МЕ ji ар| f(p)/p] il du(u* + p? — kp) (25) 
(ko— kpR)/2 nax E Vk] | 


Introducing the dimensionless variables x = p/kg, у = шір, v = kg/kg we get finally: 
5. ДЪ? Va-cvo-x 
2h?kg 1 | 
A қыл Jl ах f(x)|x] | dy(x?+y?— 1) (26) 


(»—1)2 max [(v—x), V1—x*] 


Im Рса(%) = — 


where the index Gd indicates the Goldberger method. 


$ 3. Numerical results and discussion 


In all the numerical calculations we assume the nuclear density 00 corresponding to 
the value of 1,27 f for ry = ( лоо)" (kg = 1,2 f-!) and the energy of the top of the Fermi 
sea ЛЕ = —10 MeV. This choice enables us to use the results of Verlet and Gavoret (1958). 

We will discuss the optical potential as a function of the wave number of the incoming 
nucleon. The energy E, of the incoming nucleon is connected with A, through the 
equation: 
miu n 
OM 2M 

In our computations we had to extrapolate the functions g(p) and /(р), eqs. (15). (16), 
for the energy E, > 310 MeV in the interval of about 100 MeV, since no phase shifts are 
available at these energies. In the case of the Gammel-Thaler (1957) phase shifts we encountered 
the same situation for E, < 90 MeV. However, the results are not very sensitive to these 
extrapolations because of the factors by which g and f are multiplied under the integral, 
eqs. (20), (25). 

All our numerical results are shown in Tables I, II, and III, and in Figures 1 and 2. 

Table I gives the real part of the optical potential as function of the energy Eo of the 
icoming nucleon. The G C T, GT, and SM figures were obtained with the help of eq. (20) 


E, + ЛЕ (27) 
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ТАВГЕ 1 
—Re V 
kol kp 1,5 2 3 
Eo 30 80 230 
Y (Verlet and Gavoret) 48.7 23.1 8.7 
GCT 19.3 12.0 1:2 
GT (I — 0) 2153 11.9 —2.5 
SM 22.2 14.4 1.8 
GCT 43.5 40.7 32.9 
GT all Га 39.9 31.8 16.6 
SM 32.8 27.5 15.4 
GGT 72.9 51.8 40.4 
Y+ (СТ 10 67.3 43.0 27.8 
SM 59.3 36.3 22.3 


Table І. The real part of the optical potential as function of incoming nucleon energy Ep. All figures are in MeV. 


with the Gammel-Christian-Thaler (1957), Gammel-Thaler (1957), and Signell-Marshak 
(1958) phase shifts respectively. The / values of the phase shifts considered are indicated. 
The figures Y+ ССТ (GT, SM)/> 0 are the results of adding the results of Verlet 
and Gavoret (1958), obtained with the Yamaguchi interaction, to the results obtained with 
the help of eq. (20) with the / > 0 phase shifts of Gammel, Christian, Thaler (Gammel, 


Thaler; Signell, Marshak). 


Fig. 1 shows some of the results given in Table and the experimental results taken 
from the review article Ьу Glassgold (1959). In addition, Fig. 1 shows the results obtained 
by Verlet and Gavoret with the Yamaguchi interaction 
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TABLE П 
— Im V 
КЁБ 1.5 2 3 

Е, 30 80 230 

Y (Vertlet апа Gavoret) 20.6 12.1 3.1 
GCT 22.9 9.1 1.5 

GT S 20.3 8.1 129 

5М 22.0 9.0 0.8 

GCT 9.7 4.6 12 

GT r= 8.5 3.9 I 

SM Gd 9.5 4.6 0.7 

ССТ 39.4 30.2 29.3 

GT ай [з 828 29.2 32.0 

SM 35.6 28.4 29.7 

GCT 20.8 17.0 24.2 

GT all Га 18.4 16.6 26.5 

SM Gd 20.2 16.1 24.5 

GCT eua 24.6 26.2 

Ү+ GT 12:40 30.5 24.8 28.5 
5М Gd 31.4 290-1 210 


Table II. The imaginary part of the optical potential as function of incoming nucleon energy É. All figures are 
in MeV. 
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Fig. 2. The imaginary part of the optical potential as function of the incoming nucleon energy Eo. 


252 


Table П gives the imaginary part of the optical potential as function of the energy Eo 
of the incoming nucleon. All figures were obtained analogically as in Table I. The abbre- 
viation “С d" means that eq. (25) was used instead of eq. (20). 

Fig 2 shows some of the results of Table II and the experimental results taken from 
the review article by Glassgold. The Y-curve shows the results obtained by Verlet and 
Gavoret with the Yamaguchi interaction. | 

Table ІП gives the real and imaginary part of the optical potential as function of the 
energy Еу for the case of the effective range (ER) separable interaction of eq. (12). All the 
symbols have the same meaning as in Tables I and П, except that the Yamaguchi inter- 
action is replaced by the ER -interaction. 


TABLE Ш 

ko kg | 1.5 802 3 
E, 30 | 80 230 
E R (Verlet and Gavoret) 51.5 254 13.7 
; ER+GCT (i 0) 75.7 544 45.4 
— Ке} ЕК-СТ (1-0 70.1 45.6 32.8 
ЕК+5М (0) 62.1 38.9 27.3 
WER, 23.9 16.1 6.0 
ЕВ ЕСТ > 0) 35.0 28.6 29.0 
— In (Gd) ERS GT 029 33.8 28.8 31.4 
ERED U >D 34.6 27.7 20.9 


Table ІП. Тһе real and imaginary part of the optical potential as function of incoming nucleon energy E, for 
the case of the effective range separable interaction. All figures are in MeV. 


The comparison of our results with experiments is difficult because of the uncertainity 
of the experimental result presented in the review article by Glassgold (1959). Moreover 
the experimental results give the parameters for the central part of the potential, the radial 
dependence of which has the form {1+ exp [(r — R)/a]). There might be some doubts | 
whether one can compare these experimental results with the results of our calculation for 
an infinite nuclear matter. With these reservations we compare, however, our results with 
these experimental data, especially since no other data are available at the moment. 

First of all, we notice from Table I that the GCT phase shifts give, especially at 
higher energies, a much deeper real part of the optical potential, and the comparison with 
the experimental points evidently favours the СТ and SM phase shifts. For this reason and 
since anyhow the ССТ potential was aboundoned when the nucleon-nucleon polarization 
was taken into consideration, we will not consider in our discussion the results obtained with 
the GCT phase shifts. We notice, however, that Im ZV is less senstitive to the particular 
choice of the phase shifts. On the other hand, both the GT and SM phase shifts give 
very similar results for Re V and Im V. Thus these two phase shifts turn out to be equally 
good for the description of the nucleon—nucleus scattering. 


Let us now discuss our results for the real part of the optical potential. Adding the 
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¿> 0 contribution to the results of Verlet and Gavoret increases the depth of the real part 
of the optical potential obtained by Verlet and Gavoret quite essentially, especially at higher 
energies. Similarly, the Y + GT (SM) 1 > 0 curves lay higher than the GT (SM) curves. 
The last phenomenon corresponds to the increase of the real well depth when exclusion 
principle is taken into account, as was noticed by Verlet and Gavoret (V, R > Рерріп 
their Table I). 


The experimental points seem to favour the GT (SM) curves as compared with the 
Y+ RT (SM)? > 0 curves. Even more the Y curve seems to give the best overall fit to 
the experimental points. The last fact may appear surprising because the separable nucleon- 
nucleon interaction acting in the S state only is an approximation particularly poor at high 
energies. It is known that at about 300 MeV the real part of the optical potential vanishes. 
The fact that the Y curve of Verlet and Gavoret approaches zero at high energies is simply 
the results of the assumed oversimpliefied nucleon-nucleon interaction?. Since certainly 
one cannot doubt the superiority of the calculations based on the GT (SM) phase shifts, 
one has to look for other reasons for the vanishing of Re V at about 300 MeV. However, 
atthe moment we do not know these reasons. The same difficulty was encountered by 
Lippmann and Moravesik (1960) in their optical model calculations based on the recent 
modifications of the analysis of nucleon-nucleon scattering (Cziffra at al. 1959)2. 


As was mentioned in the introduction, our calculations are not applicable at low energy. 
Still, let us make one comment on the behaviour of Re V calculated according to eq. (19). 
For ky = kp the u-integration gives zero, Since, on the other hand, g(x) is regular at 
x = 0 (constant S-wave term at zero energy) we get Ке V = 0 for ky = kp, in disagree- 
ment with current ideas. However, Re V calculated with the t matrix obtained by solving 
eq. (2) with the exclusion principle taken into account, does not vanish at kọ = kg as 
shown by the results of Verlet and Gavoret. This indicates that the difficulty at ky = kp 
noticed by Lippman and Moravesik (1960) may be removed by using the ¢ matrix rather 
than the tọ matrix. 

Let us now discuss our results for the imaginary part of the optical potential. Similarly 
as in the case of Re V, adding the / > 0 contributions to the results of Verlet and Gavoret 
increases the depht of the imaginary part of the optical potential calculated by Verlet and 
Gavoret. At first sight it seems surprising that the curves Y + СТ (SM) ¿> 0 (Gd) lay 
higher than the curves G T (S M) (all J's, Gd). To explain it let us notice that the difference 
between the two groups of curves is equal to the difference in the corresponding S-wave 
contributions. In the case of the Y+ GT (SM) / > 0 (Gd) curves the 5 wave contri- 
bution is calculated by solving exactly eq. (2) for the ¿ matrix (although for the simpli- 
fied Yamaguchi potential), whereas in the case of the G T (S M) (ай /s, Gd) curves the 5 
wave contribution is calculated with the help of the Goldberger method. However, it is 


2 The same remark concerns the results obtained by Riesenfeld and Watson (1956) with help of the Fesh- 
bach-Lomon (1956) phase shifts. | 

3 The deviations from the Fermi model of the nucleus, i. e. the nucleon-nucleon correlations in the target 
nucleus give corrections diminishing the depth of the real part of the optical potential. These corrections, 


however, seem to be much too small to explain the present discrepancy (Dabrowski and Sawicki 1959). 
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known (Verlet and Gavoret 1958) that the Goldberger method badly overestimates the role 
of the Pauli principle in reducing Im V at low and intermediate energies. This is then the 
reason for our behaviour of the Y-- GT (SM) / — 0 (Gd) curves as.compared with the 
G T (SM) (all /5, Gd) curves. 

According to the above remarks the good agreement with experiment shown by the 
GT (SM) (all l's, Gd) curves seems to be accidental. One should rather look for factors 
which could reduce the Y+ G T (SM) ¿> 0 (Gd) results. 

One may notice in Table II that the modification of the / > 0 contribution introduced 
by the exclusion principle is not essential, especially since the Goldberger method over- 
estimates the effect. This supports our “Y+ GT (SM) / > 0" method of calculating the 
optical potential. 

The curves in Fig. 2 are drown even at low energies, where our theory is not aplicable 
Of course, Im V = 0 at ky = kg because of the presence of the Q operator in the Goldberger 
method in eq. (25). However, even if we used eq. (19) rather than eq. (25), thus disregarding 
the exclusion principle, we still would get Im V = 0 at ky = kg for exactly the same reason 
for which Re V = 0 according to eq. (19) at ky = Ёк, as it was mentioned above. The 
vanishing of Im V at k, = kg certainly agrees with current notions. 

We cannot compare our results for Im V with the results of Gomes (1959). By using the 
independent pair wave function calculated by Gomes at al. (1958), he restricts himself to 
low energy. For reasons, however, mentioned in the introduction we do not feel we are 
able to apply our expressions for Im V for low energy. In contradistinction to the procedure 
of Gomes, by using the Yamaguchi interaction one does not neglect the singlet-triplet 
difference in the S state nucleon-nucleon interaction. 

Table ПІ repeats the results of Table I and П in the case of the effective range separ- 
able nucleon-nucleon potential. The results in Table III do not differ esentially from our 
previous results. 

Тһе authors express their gratitude to dr Gomes and to dr Moravesik for communica- 
tion of their results prior to publication. 
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ENERGY LEVELS IN :%Ег NUCLEUS 
By E. Bozex, Н. NIEWODNICZANSKI, S. Ocaza, S. SzymczyK AND Yu. V. NorsEEv* 
Institute of Nuclear Physics, Cracow 
(Received October 7, 1960) 


The level scheme of 186Ег has been investigated from the decay of 188Тш (7.7 h) іп equi- 
librium with 186ҮЬ (58h). Gamma-gamma coincidence experiments were performed and a new 
level scheme for 198Ег is proposed with levels at 80, 264, 788, 862, 957, 1536, 2140 and 2354 keV 
with corresponding spin and parity assignments. 


1. Introduction 


The purpose of this work was to obtain more informations concerning the level scheme 
of 166Er resulting from the decay of !99Yb. 

The !99Yb (58 h) nucleus decays by electron capture to 166Tm (7.7 h), which decays 
by electron capture and a small amount of В+ to Er. 

In the course of the last few years the energies of gamma radiation and of conversion 
electrons in the decay of 166Yb and 16Tm were measured by different authors [1—6]. 

An incomplete level sheme of 16 Ет from the decay of 166Tm was proposed by Gromov 
et al. [7] and Jacob et al. [2]. Their proposals were based on measurements of energies and 
intensities of transitions. Boskma and De Waard [6] on the basis of conversion electron- 
gamma and gamma-gamma coincidence experiments proposed a level scheme for the 
most intensive transitions іп 19Ет. 

In the present paper some gamma-gamma coincidence experiments in the decay of 
166ү} are described and a more complete level scheme of 1%Er is proposed. 


2. Sources and. Apparatus 


The isotopes of Yb were obtained in the reaction of spallation by 3-4 hours irradiation 
of Ta targets with 660 MeV protons from the synchrocyclotron of the Joint Institute for 
Nuclear Research in Dubna (USSR). The chemical operations — the isolation of the rare 
earths group from the Ta target and the chromatographic separation of elements — were 
performed in Dubna 11-19 hours after the irradiation. 


* Joint Institute for Nuclear Research, Dubna. 
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The measurements in Cracow started 1.5 to 3 days after each irradiation. In the experi- 
ments reported in this paper the samples of four Ta targets were used. The irradiations were: 
performed in about 2 months intervals. In the first two runs the measurements of the 
y-radiation were made at different times after the irradiation for the decay analysis 
of the entire spectrum of the isotope mixture of ytterbium. In the next two runs the gamma- 
gamma coincidence measurements were performed. 

The single count y-spectrum was recorded with a scintillation spectrometer using: 
a 13" x 1" NaJ(TI) crystal, an EMI 6097 F photomultiplier and а 100-channel pulse: 
height analyser. The resolving power of this spectrometer for the 661 keV 187Cs line was: 
8.290. 

The gamma-gamma coincidence spectra were recorded on the multichannel pulse: 
height analyzer, gating the spectrum from one counter with coincidence pulses from: 
a fast-slow coincidence circuit. The block scheme of the apparatus is shown in Fig. 1.. 


` 


^a (7t ) 


/77-сдате/ 
Analyzer 


Fig. 1. Block diagram of electronics. 


NaJ (TI) crystals 12" x 2" and 2" x 2" and ЕМІ 6097 F photo-multipliers were used. The: 
resolving powers of both spectrometers for the 137Cs line were 9%. The resolution time: 
of the coincidence cireut 2r = 5.5 x 10-8 sec. 


3. Results of Measurements 


One could expect the investigated sources to contain radiation of the following neutron: 
deficite isotopes: 166ү} in equilibrium with Tm, 167YD with its daughter nucleus !9"Tm, 
and 1€? yp. 

The single count y-spectrum recorded on the 100-channel amplitude analyzer 
is shown in Fig. 2a and 2b. 'The decay analysis of the photopeaks shows, that they decay with. 
a half-period of about 57 hours, i. e. they belong to the spectrum of 16ҮЪ, No activity with. 
a half-period of 9.6 days was found, which could belong to the activity of 19"Tm, proceeding | 
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from the decay of 167ҮЪ (18.5 min) and 167ҮЬ* (74 min). This сап be explained by the time 
interval between the end of irradiation of Ta and the moment of separation of the ytter- 
bium sample from the other rare earth elements. Few days after the measurements started 


the decay of the 155, 180, 250 and 300 keV lines with a half-period of 30.6 days could be 
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Fig. 2a. Low energy part of the y-spectrum of Yb. 
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Fig. 2b. High energy part of the y-spectrum of Yb. 
observed. This is the decay of 169 Yb. From this decay analysis it was established, that the 


coincidence experiments for 166Yb may be carried out in the first week of measurements, 
when the 166Yb activity is the strongest compared to the other present in the sample. The 
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comparison of the energies of y-radiation of 198 Ет obtained in our measurements withl 
those published by other authors is shown in Table I. 

The existence of a transition with the energy 880 keV reported up to now only by 
Brabec et al. [4] is confirmed by our coincidence measurements. 


TABLE I 


Summary of observed y-rays (in keV) 


Magnetic spectrometers Scintillation spectrometers 
Gromov Mihelich Boskma Boskma Baranovskii Cat adum 
[7] [1] [6] [6] [5] 
80.0 807 | 80.2 80 80 80 
153.4 154.6 — 153 
183.9 184.7 183 180 184 
193.6 194.8 192 194 185 
214.1 215.4 = 215 
282 
299 
347 
405 
433 
459 464 464 450 460*) 
522 522 
596 596 600 
598 
674 
693 710 710 690 710 
707 
749 
760 
782 790 790 780 790 
789 | 
(807) 
877 880*) 
(1158) 
1181 1180 1180 1180 
(1209) 
1277 1270 1280 1270 
(1309) 
1380 
(1423) 1380 (1380) 
! Ue 1550 
2060 2000 3060 es 
2088 060 2070 


۹ : к 
) Observed in the coincidence spectrum 
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We have not found neither the lines 153 and 522 keV measured with scintillation 
spectrometer by Baranovski and Pokrovskii [5], nor the line 1550 keV reported by Boskma 
апа De Waard [6]. š 


The widths and shapes of photopeaks with energies 185, 710, 790 keV show, that 
[һезе peaks contain more than one transition. 


TABLE II 


Coincidence relationships of lines in !99Yb decay 


80 185 | 460 | 600 | 710 | 790 | 880 | 1180 | 1270 | 1880 | 2070 
185 | ы | + + = = + - 
ES | - 4 rM e (+) | | E 
790 - | | + = = d + = 5 
1180 4 1 = = 5 : 
1270 + | 24 = | + + | — = -4 =. - 
2070 4 + -=% | - - - - 

| 


| The results of our coincidence measurements are listed in Table II. The lines 185, 
710, 790, 1180, 1270 and 2070 keV were discriminated by a single channel pulse height 
analyzer. The measurements for each line were performed applying different amplifications 
of the amplifiers, different channel widths and different positions of the channel on the peak. 
[In this manner both the energy intervals of the coincidence spectrum recorded by the multi- 
channel analyzer and the energy intervals selected by the single channel analyzer were 
changed. 

The following coincidences stated in the present work are in agreement with the results 
of Boskma and De Waard [6]: the 185 keV line with photopeaks 80, 710, 790, 1180 keV; 
the 710 keV line with 80, 185, 1180 keV peaks; the 790 keV line with 80 and 1270 keV 
peaks; the 1180 keV line with 80, 185, 710 keV peaks and the 1270 keV line with 80 and 
790 keV peaks. 

Besides this confirmation of the results of Boskma and De Waard [6] we obtained in 
coincidence with the 185 keV line the photopeaks 185, 460 and 1880 keV (Fig. 3a, b), but 
sould not confirm the coincidence of the 185 keV line with the 2070 keV peak reported 
by Boskma and De Waard. Further we obtained the following coincidences: between the 
710 keV line and 460, 710, 790 keV peaks (Fig. 4); between the 790 keV line апа 185, 
160, 710 keV peaks (Fig. 4); between the 1180 keV and the 880 keV peak (Fig. 5); between 
һе 1270 keV line and 185, 600 keV peaks (Fig. 5) and between the 2070 keV line and 80, 
185 keV photopeaks. None of these coincidence relationships was reported by Boskma and 
De Waard. 


4. Level Scheme 


On the basis of our experiments and the energy values of transitions given by Gromov 
t al. [7] (see Table I) we are suggesting a level scheme for 166Er from the decay of 169Tm. 


[his level scheme is shown in Fig. б. 
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Fig. 3a. ү-ү coincidences of the 460, 603 and Fig. 3b. y-y coincidences of the 1880 keV line with the 
710 keV lines with the 185 keV line in 19Ет, 185 keV line in 166Fr. 
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Fig. 4. ү-ү coincidences with the 710 and 790 keV lines іп 166Er 
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Fig. 6. The level scheme of !66Er. 


264 


The two levels at 80 and 264 keV were reported by Boskma and De Waard [6] ап 
are also known from the 166Ho decay. They are interpreted as rotational states in agreemer 
with the #2 character of 80 and 183,9 keV transitions. Their assignments are (0,2 +) аш 
(0,4 +) respectively. 

The three levels at 788, 862 and 957 keV сап be interpreted as belonging to a y vibra 
tional state and to two lowest rotational states built on it. Their assignments would be (2,2 + 
(2,3 +) and (2,4 +) respectively. Boskma and De Waard have interpreted the level 862 ke" 
as a vibrational state with (2,2 +) and the level at 957 keV as a rotational state (2,3 +). The: 
did not mention the 788 keV level. In the level scheme proposed by Boskma and De Waar: 
there is an inconsistency, even regarded so by them. According to their scheme two transi 
tions with about the same energy of 870 keV must occur, the first between the 862 keV 
level and the ground state, the second one between the 957 and 80 keV levels. The tota 
intensity of both transitions would be sufficiently large to give an observable photopeali 
in the single count gamma spectrum. In our proposal (suggested earlier by Jacob et al. [2: 
and Gromov et al. [7], but not supported by any coincidence measurements) this inconsis: 
tency is removed. In Table IIT there are listed the branching ratios of transitions betweer 


TABLE III 


The branching rations of trasitions from levels with K — 2 to lewels with K — 0 in 160 Ey 


Branching ratios 
Levels and š 
D шар 2 theoretical 
transitions experimental [7] 
Boskma, De Waard [6] Our proposal 
788 keV level 
7 A. I, (2—0 
y, 789 keV — = 0.63 = | A Bae Ky 
Ya 707 keV 1 | L, (2 > 2) 
862 keV level i ne Di (2 > 2) 635 T. (3--2 95 
y, 782 keV Laka L(2—4) 7 L, (3 — 4) 5 
Ya 598 keV 3 Š 2 
th oe ,ر1‎ (2 2) 2, I,,, (8 — 2) < 
Lm 1з (2 ج‎ 0) 1,»(8 — 0) 
957 keV level 
y, 817 keV 1.2 0.66 D,-2) _ m 4-2 _ 
Уз 693 keV 1, I, (8 — 4) 1, (4-> 4) 


the levels with K = 2 and levels with K = 0 calculated on the ground of Alaga's theory for 
the level scheme proposed by Boskma and De Waard and for the level scheme suggested 
by us. For comparison the table contains also the experimental results given in the paper 
of Gromov et al. [7]. From this table one can clearly see, that the level scheme proposed 
by us is in better agreement with the experimental branching ratios. 

The coincidence relationship between the photopeaks at 710 and 190 keV can be 
explained by the existence of a level at 1536 keV. Because of the El character of the 674 and 
749 keV transitions determined by the conversion coefficients measurements [7], the level 
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at 1536 keV may have ап assignment (2,2—). However, the experimentally determined 
branching ratio of these transitions is in disagreement with the theoretically expect- 
ed one. 

The level at 2140 keV is also reported by Boskma and De Waard, who assigned it 
(2—) on the ground of the log ft calculation of the electron capture transition from the ground 
state of 199Tm to this state and taking into account the measured banching ratio of 1180 an 
1270 keV transitions. On the other side, Gromov et al. [7] reported, that both, the 1181 and 
the 1277 keV transitions are of E2 or possibly (M1 + E2) type, and hence the level at 2140 keV 
would have the positive parity. As no transition with the energy 1352 keV, between the levels 
at 2140 and 788 keV, was found, the spin of the 2140 keV level would be at least 5. This 
is, however, rather improbable due to the small value of log ft of the electron capture transition 
to this level [6] together with the spin 2 or 3 and negative parity) of the ground state of 
166[m, The latter values of spin and parity are expected on the ground of Nilson's 
caleulations of single particle levels in deformed nuclei and the comparison evidence of 
spins and parities of the neighbouring even-odd nuclei with the same number of protons 
or neutrons. 

Тһе level at 2354 keV may be the same to which Boskma and De Waard ascribed the 
energy 2260 keV basing on the coincidence relationship 180-2000 keV. The energy 2354 keV 
ascribed to this level allows to introduce in the level scheme the 214.1 keV transition. If 
we put this transition between the levels 2354 and 2140 keV, then the coincidence relation- 
ships between the 185 keV photopeak and each of the discriminated energies 185, 790 and 
1270 keV can be explained. The Е1 character of the 214.1 keV transition [7] in connection 
with the positive parity of the 2140 keV level requires the negative parity for the 2354 keV 
level in agreement with the suggestion of Boskma and De Waard. In our measurements 
we have not obtained the 2070 keV photopeak in the coincidence spectrum with the line 
185 keV, but otherwise in the coincidence spectrum with the line 2070 keV we have obtained 
the 80 and 185 keV lines. This fact can possibly be explained if we accept, that only the 
transition 2088 keV is in coincidence with the 183.9 keV line and that the transition 2060 keV 
(observed in the conversion electron spectrum [4]) is that between the 2140 and 80 keV 
levels. 

The 460 keV line which was observed in our measurements in the coinci- 
dence spectra with 185, 710 and 790 keV lines is not introduced into the proposed level 
scheme. This can not be a transition between the levels proposed in the present 
paper. In order to explain the appearance of this line a new level is required. However, 
for the proper introduction of this level into the level scheme our informations are not 
sufficient. 

Our proposal of the level scheme for 199Er needs some more confirmation. To 
get this, we are starting now the y-y directional correlation. experiments in order to 
obtain definite values of spins for the 788, 862, 957 and 2140 keV levels. 

The authors are indebted to Mr 5. Chojnacki and Mr I. A. Yutlandov for the 
work done with the preparation of sources. Thanks are also due to Mr K. Ma- 
linowski and Mr T. Walczak for their valuable assistance in the course of meas- 


urements. 
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Fig. 1. The influence of 
density fluctuations оп 
ultrasonic waves in the 
solution of aniline in cyclo- 
hexane near its critical 
point for several temper- 
atures: 
азд (бе (02-21 CE 
с 4 30199 (rs 
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Fig. 2. The ultrasonic field and density fluctuations in the solution of aniline in cyclohexane at 30.9°С (three- 


fold magnified). 
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INFLUENCE OF DENSITY FLUCTUATIONS ON ULTRASONIC WAVES IN A BINARY 
LIQUID MIXTURE NEAR ITS CRITICAL SOLUTION POINT 


Bv ANTONI ŜLIWIŃSKI 
Institute of Acoustic and Vibration Theory, Mickiewicz University, Poznañ 


(Received December 8, 1960) 


Brown and Richardson (1959) measured the scattering of the ultrasonic radiation 
in aniline-cyclohexane mixture near the critical solution point. They also observed optically 
the fluctuation clusters by means of the schlieren method. 

It is interesting to observe during a longer time and at various temperature the influence 
of the density fluctuations on ultrasonic waves passing through the mixture near its critical 
point. The approaching to the critical solution temperature when the temperature is lowered 
step by step is shown in Fig. 1; photographs were obtained by the schlieren technique. 
There are thirteen momentary pictures made at constant temperatures in time intervals 
of 1/12 sec between each other, arranged in seven columns. In photographs a, b, c, e, g both 
clusters and fronts of the standing ultrasonic wave appear simultaneoulsy. In photographs d, f, 
the ultrasonic field was switched off and these pictures are given only to be compared with 
pictures c, е. In picture g the ultrasonic field was switched on but the fluctuations became 
so strong that the wave fronts could not be seen. All the pictures in Fig. 1 are reduced in 
size. Fig. 2 shows one photograph in threefold magnification. 

As it may be seen from Fig. 2 the wave fronts were broken off at many places. When 
the temperature decreased approaching the critical point these breaches became frequent 
and simultaneously the sizes of the clustres decreased. Sometimes one сап observe a “brick 
wall structure" mentioned by Ernst (1947). 

Тһе size of the clusters at a given temperature may be determined from the photo- 
graphs directly by comparing it with the length of the ultrasonic wave. The number of 
breaches of the wave fronts in a given area of ultrasonic field, e.g. in 4 cm? may be considered 
as the measure of the crumbling of the fluctuations (sizes of clusters). In this way one can 
obtain the crumbling of fluctuations as a function of temperature, counting the number 
of breaches of wave fronts in the various photographs at various temperatures (see Fig. 3). 
As it may be seen there is a large dispersion between the series calculated for various repeti- 
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tions of the experiments. Exactly at the critical point (30.7°С) the process of crumbling 
of clusters was so far advanced that the medium was looking like a fog. 

The observations were made in a 33.3 per cent solution of aniline in cyclohexane. 
The critical solution temperature of the mixture was 30.7°C. The frequency of ultrasonic 
was 1 MHz, the temperature was kept constant within 0.05*C. 


300 


250 


ho 
© 
© 


Number of breaches of wave fronts in 4 cm? 


E 
сч 
© 


100 


50 


t, 2307 31 315 32 32.5 C 
temperature 


Fig. 3. The crumbling of fluctuations in the solution of aniline in 


tion of temperature. 


cyclohexane near its critical point as a func- 
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ON THE EFFECT OF THE AGEING TIME OF THE DEVELOPPER ON NUCLEAR 
EMULSION BACKGROUND 
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(Received May 25, 1960; revised manuscript received January 17, 1961) 


The present paper deals with the effect of the time of ageing of the developper on the 
background of the nuclear emulsion. As time proceeds, the developper will be undergoing 
oxidation, and so its effect will gradually become weaker. Hence, in developping a nuclear 
plate with an “old” developper, the background should be expected to be lesser than in one 
developped with a “fresh” developper. 

Nuclear “Agfa” K-2 plates 400ш in thickness were used. These were saturated with 
lithium citrate and then irradiated with slow neutrons. This yielded 


SLi 4- in He 4- 7T. 


An equal number of plates were developped by the temperature method (amidol developper) 
and by the method of decreased PH (Lohman 1956). In proceeding by the latter method, 
the composition of the developper and the stage 1 to stage 2 time ratio were varied. On the 
plates developped by either method, the paths (T --He) are well visible. The background, 
however, changes according to the time of ageing of the developper. In order to obtain 
a quantitative relationship, countings were made of the number of grains of the background 
per unit area of the emulsion, for different periods of ageing of the developper. From the 
data thus obtained, the value of the background T was plotted versus the time t of ageing 
of the developper. These are shown in Fig. 1 for plates developped by the temperature 
method, and in Fig. 2 for those developped by the method *of lowered PH." 

From Figs. 1 and 2, the background in plates developped with the “fresh” developper 
(aged for 0—1.5 days) is denser than in those developped with “old” developper (time of 
ageing exceeding 1.5 days). 

The experiments lead to the following conclusions: 1. One should avoid developping 
nuclear plates with either too fresh or too old a developper. The authors propose the time 
interval of 1.5 to 3 days for the method “of lowered РН” and over 2 days for the tempera- 
ture method. These intervals provide for a constant background value. 2. The background 
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s considerably lesser in plates developped by the method “of lowered PH? (approximately 
ne tenth). 
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